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A variety of observations suggest that magnetic fields are present in all galaxies and galaxy clusters. 
These fields are characterized by a modest strength (10~^ — 10~^ G) and huge spatial scale (< 
IMpc). It is generally assumed that magnetic fields in spiral galaxies arise from the combined 
action of differential rotation and helical turbulence, a process known as the «a;-dynamo. However 
fundamental questions concerning the nature of the dynamo as well as the origin of the seed fields 
necessary to prime it remain unclear. Moreover, the standard atJ-dynamo does not explain the 
existence of magnetic fields in elliptical galaxies and clusters. The author summarizes what is 
known observationally about magnetic fields in galaxies, clusters, superclusters, and beyond. He 
then reviews the standard dynamo paradigm, the challenges that have been leveled against it, and 
several alternative scenarios. He concludes with a discussion of astrophysical and early Universe 
candidates for seed fields. 
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I. INTRODUCTION 

The origin of galactic and extragalactic magnetic fields is one of the most fascinating and challenging problems 
in modern astrophysics. Magnetic fields are detected in galaxies of all types and in galaxy clusters whenever the 
appropriate observations are made. In addition there is mounting evidence that they exist in galaxies at cosmological 
redshifts. It is generally assumed that the large-scale magnetic fields observed in disk galaxies are amplified and 
maintained by an aw-dynamo wherein new field is regenerated continuously by the combined action of differential 
rotation and helical turbulence. By contrast, the magnetic fields in non-rotating or slowly rotating systems such as 
elliptical galaxies and clusters appear to have a characteristic coherence scale much smaller than the size of the system 
itself. These fields may be generated by a local, turbulent dynamo where, in the absence of rapid rotation, the field 
does not organize on large scales. 

In and of itself, the dynamo paradigm must be considered incomplete since it does not explain the origin of the 
initial fields that act as seeds for subsequent dynamo action. Moreover, the timescale for field amplification in the 
standard aw-dynamo may be too long to explain the fields observed in very young galaxies. 

It is doubtful that magnetic fields have ever played a primary role in shaping the large-scale properties of galaxies 
and clusters. In present-day spirals, for example, the energy in the magnetic field is small as compared to the rotation 
energy in the disk. To be sure, magnetic fields are an important component of the interstellar medium (ISM) having an 
energy density that is comparable to the energy density in cosmic rays and in the turbulent motions of the interstellar 
gas. In addition, magnetic fields can remove angular momentum from protostellar clouds allowing star formation to 
proceed. Thus, magnetic fields can play a supporting role in the formation and evolution of galaxies and clusters but 
are probably not essential to our understanding of large-scale structure in the Universe. 

The converse is not true: An understanding of structure formation is paramount to the problem of galactic and 
extragalactic magnetic fields. Magnetic fields can be created in active galactic nuclei (AGN), in the first generation 
of stars, in the shocks that arise during the collapse of protogalaxies, and in the early Universe. In each case, one 
must understand how the fields evolve during the epoch of structure formation to see if they are suitable as seeds 
for dynamo action. For example, magnetic fields will be amplified during structure formation by the stretching and 
compression of field lines that occur during the gravitational collapse of protogalactic gas clouds. In spiral galaxies, 
for example, these processes occur prior to disk formation and can amplify a primordial seed field by several orders 
of magnitude. 

In principle, one should be able to follow the evolution of magnetic fields from their creation as seed fields through 
to the dynamo phase characteristic of mature galaxies. Until recently, theories of structure formation did not possess 
the sophistication necessary for such a program. Rather, it had been common practice to treat dynamo action and 
the creation of seed fields as distinct aspects of a single problem. Recent advances in observational and theoretical 
cosmology have greatly improved our understanding of structure formation. Ultra-deep observations, for example, 
have provided snapshots of disk galaxies in an embryonic state while numerical simulations have enabled researchers 
to follow an individual galaxy from linear perturbation to a fully-evolved disk-halo system. With these advances, a 
more complete understanding of astrophysical magnetic fields may soon be possible. 

This review brings together observational and theoretical results from the study of galactic and extragalactic 
magnetic fields, the pieces of a puzzle, if you like, which, once fully assembled, will provide a coherent picture of 
cosmic magnetic fields. An outline of the review is as follows: In Section II we summarize useful results from 
magnetohydrodynamics and cosmology. Observations of galactic and extragalactic magnetic fields are described in 
Section III. We begin with a review of four common methods used to detect magnetic fields; syncrotron emission, 
Faraday rotation, Zeeman splitting, and optical polarization of starlight (Section III. A). The magnetic fields in spiral 
galaxies, ellipticals, and galaxy clusters are reviewed in Sections III.B-III.D while observations of magnetic fields 
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in objects at cosmological redshifts are described in Section III.E. The latter are essential to our understanding of 
the origin of galactic fields since they constrain the time available for dynamo action. Section III concludes with a 
discussion of observational limits on the properties of cosmological magnetic fields. 

Magnetic dynamos are discussed in Section IV. We first review the primordial field hypothesis wherein large scale 
magnetic fields, created in an epoch prior to galaxy formation, are swept up by the material that forms the galaxy 
and amplified by differential rotation. The model has serious flaws but is nevertheless instructive for the discussion 
that follows. Mean-field dynamo theory is reviewed in Section IV. B. The equations for a disk dynamo are presented 
in Section IV. C and a simple estimate for the amplification rate in galaxies is given in Section IV. D. 

The standard mean-field treatment fails to take into account backreaction of small-scale magnetic fields on the 
turbulent motions of the fluid. Backreaction is a potentially fatal problem for the dynamo hypothesis for if magnetic 
fields inhibit turbulence, the dynamo will shut down. These issues are discussed in Section IV. E. 

Galactic magnetic fields, like galaxies themselves, display a remarkable variety of structure and thus an under- 
standing of galactic dynamos has required full three-dimensional simulations. Techniques for performing numerical 
simulations are reviewed in Section IV.F and their application to the problem of diversity in galactic magnetic fields is 
discussed in Section IV. G. In Section IV. H we turn to alternatives to the acj-dynamo. These models were constructed 
to address various difhculties with the standard scenario. Section IV ends with a brief discussion of the generation of 
magnetic fields in elliptical galaxies and galaxy clusters. 

The question of seed fields has prompted a diverse and imaginative array of proposals. The requirements for seed 
fields are derived in Section V.A. Section V.B describes astrophysical candidates for seed fields while more speculative 
mechanisms that operate in the exotic environment of the early Universe are discussed in Section V.C. 

The literature on galactic and extragalactic magnetic fields is extensive. Reviews include the excellent text by 
Ruzmaikin, Sokoloff, & Shukurov (1988a) as well as articles by Rees (1987), Kronberg (1994), and Zweibel & Heiles 
(1997). The reader interested in magnetohydrodynamics and dynamo theory is referred to the classic texts by Moffatt 
(1978), Parker (1979), and Krause & Radler (1980) as well as "The Almighty Chance" by Zel'dovich, Ruzmaikin, 
& Sokoloff (1990). A survey of observational results from the Galaxy to cosmological scales can be found in Vallee 
(1997). The structure of galactic magnetic fields and galactic dynamo models are discussed in Sofue, Fujimoto & 
Wielebinski (1986), Krause & Wielebinski (1991), Beck et al. (1996), and Beck (2000) as well as the review articles 
and texts cited above. 



II. PRELIMINARIES 

A. Magnetohydrodynamics and Plasma Physics 

Magnetohydrodynamics (MHD) and plasma physics describe the interaction between electromagnetic fields and 
conducting fluids (see, for example, Jackson 1975; Moffatt 1978; Parker 1979; Freidberg 1987; Sturrock 1994). MHD 
is an approximation that holds when charge separation effects are negligible. Matter is described as a single conducting 
fluid characterized by a density field /9(x, t), velocity field V(x, t), pressure p(x, t), and current density J(x, t). The 
simple form of Ohm's law is valid while the displacement current in Ampere's Law is ignored. In Gaussian units, the 
relevant Maxwell equations take the form 



V-B = 



(1) 




(2) 




(3) 



and Ohm's law is given by 



J' 



crE' 



(4) 



where a is the conductivity and "primed" quantities refer to the rest frame of the fluid. Most astrophysical fluids are 
electrically neutral and nonrelativistic so that J' = J and E' = E + (V x B) /c. Eq. (^) becomes 
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J = o-fE+ — — j (5) 
which, when combined with Eqs. (||) and (^), yields the ideal MHD equation: 

— = V X (V X B) + TyV^B . (6) 

In deriving this equation, the molecular diffusion coefficient, rj = c^/Aira, is assumed to be constant in space. 
In the limit of infinite conductivity, magnetic diffusion is ignored and the MHD equation becomes 

— = V X (V X B) (7) 

or equivalently 

— = (B.V)V-B(V.V) (8) 

where d/dt — d / dt + \ ■ 'V is the convective derivative. The interpretation of this equation is that the flux through 
any loop moving with the fluid is constant (see, for example, Jackson 1975; Moffatt 1978; Parker 1979), i.e., magnetic 
field lines are frozen into the fluid. Using index notation we have 

dt dxj * dxj 

where a sum over repeated indices is implied. This equation, together with the continuity equation 

f = -,V.V, ,10, 

gives 

dB, 2B,dp 

where ct^ = djVi — \5ijdkVk (see, for example, Gnedin, Ferrara, & Zweibel 2000). 

The appearance of convective derivatives in Eq. (|l|) suggests a Lagrangian description in which the field strength 
and fluid density are calculated along the orbits of the fluid elements. The first term on the right-hand side of Eq. ( pT| ) 
describes the adiabatic compression or expansion of magnetic field that occurs when V • V 7^ 0. Consider, for example, 
a region of uniform density p and volume V that is undergoing homogeneous collapse or expansion so that Uij = 
and V • V = C where C = Cit) is a function of time but not position. Eq. ( pT|) implies that B oc p^/"^ cx V^^^^ . 
Thus, magnetic fields in a system that is undergoing gravitational collapse are amplified while cosmological fields in 
an expanding universe are diluted. 

The second term in Eq. ( pi] ) describes the stretching of magnetic field lines that occurs in flows with shear and 
vorticity. As an illustrative example, consider an initial magnetic field B — Bo's, subject to a velocity field with 
dVy/dx — constant. Over a time t, B developes a component in the y-direction and its strength increases by a factor 



1/2 

Combining Eq. (^) and (|l^) yields the following alternative form for the MHD equation: 



l + {tdVy/dxf 
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The formal solution of this equation is 



where ^ is the Lagrangian coordinate for the fluid: 



X^{t)^C^+ f V,{s)ds . (14) 

Jo 

It follows that if a "material curve" coincides with a magnetic field line at some initial time then, in the limit rj = 0, 
it will coincide with the same field line for all subsequent times. Thus, the evolution of a magnetic field line can be 
determined by following the motion of a material curve (in practice, traced out by test particles) as it is carried along 
by the fluid. 

The equation of motion for the fluid is given by 



— + (V-V)V = — Vp- V$ + — (J X B) + i/V2V (15) 
ot p cp 

where h' is the viscosity coefficient and $ the gravitational potential. In many situations, the fields are weak and the 
Lorentz term in Eq. ( p^ ) can be ignored. This is the kinematic regime. In the limit that the pressure term is also 
negligible, the vorticity ^ = V x V obeys an equation that is similar, in form, to Eq. (H): 



dC 

^ = Vx (VxO + i^V^C . (16) 



Moreover, if viscosity is negligible, then ^ satisfies the Cauchy equation (Moffatt 1978): 



Q (x, t) ^ Q (I, 0) dx, 

p{^,t) p(e 0) ae, ■ ^ ' 

However, Eq. (^^ is not a solution of the vorticity equation so much as a restatement of Eq. ( pj| ) since dx/dS, is 
determined from the velocity field which, in turn, depends on x. By contrast, in the kinematic regime and in the 
absence of magnetic diffusion, Eq. (|l^ provides an explicit solution of Eq. (pT|). 

The magnetic energy density associated with a field of strength B is es = B'^/Stt. For reference, we note that the 
energy density of a 1 G field is ~ 0.040 erg cm^'^. A magnetic field that is in equipartition with a fiuid of density p 

1/2 

and rms velocity v has a field strength B ~ B^q = (47rpi'^) . In a fiuid in which magnetic and kinetic energies are 

1 /2 

comparable, hydromagnetic waves propagate at speeds close to the so-called Alfven speed, va = {B'^/iirp) 

It is often useful to isolate the contribution to the magnetic field associated with a particular length scale L. 
Following Rees & Reinhardt (1972) we write 



B^\ _ f B{LfdL ^^g^ 



Stt / Stt L 

where (B^/Stt) is the magnetic field energy density averaged over some large volume V. B{L) is roughly the 

1/2 

component of the field with characteristic scale between L and 2L. Formally, B(i) — (fc'^/27r^V) where 
Bfc = / cfixexp (ik • x)B(x) is the Fourier component of B associated with the wavenumber k = 2tt/L. 

In the MHD limit, magnetic fields are distorted and amplified (or diluted) but no net fiux is created. A corollary of 
this statement is that if at any time B is zero everywhere, it must be zero at all times. This conclusion follows directly 
from the assumption that charge separation effects are negligible. When this assumption breaks down, currents driven 
by non-electromagnetic forces can create magnetic fields even if B is initially zero. 
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B. Cosmology 

Occasionally, we will make reference to specific cosmological models. A common assumption of these models is that 
on large scales, the Universe is approximately homogeneous and isotropic. Spacetime can then be described by the 
Robertson- Walker metric: 



ds^ :=^c^df ~a^{t)dr^ (19) 

where a(t) is the scale factor and dr is the three-dimensional line element which encodes the spatial curvature of the 
model (flat, open, or closed). For convenience, we set a(fo) — 1 where is the present age of the Universe. The 
evolution of a is described by the Friedmann equation (see, for example, Kolb & Turner 1990) 

SttG , , k 

= ^5- er + em + CA 20 

where e^, Cm, and eA are the energy densities in relativistic particles, nonrelativistic particles, and vacuum energy 
respectively, k = 0, ±1 parametrizes the spatial curvature, and H is the Hubble parameter. Eq. ( po| ) can be recast as 



1 /2 

Ida rr f , , ^ , {I - ~ ^rn - ^a)\ ^ 
--T7=-Wo — H 3-+i2AH 5 (21) 

a dt \ a* a"^ a- J 

where Hq = H{t()) is the Hubble constant and is the present-day energy density in units of the critical density 
Cc = 3H'^/8ttG, i.e., fl^ = Em/ec, etc.. 

Recent measurements of the angular anisotropy spectrum of the CMB indicate that the Universe is spatially flat 
or very nearly so (Balbi et al. 2000; Melchiorvi et al. 2000; Pryke et al. 2001). If these results are combined with 
dynamical estimates of the density of clustering matter (i.e., dark matter plus baryonic matter) and with data on 
Type la supernova, a picture emerges of a universe with zero spatial curvature, — 0.15 — 0.4, and = 1 — il™ 
(see, for example, Bahcall et al. 1999). In addition, the Hubble constant has now been determined to an accuracy 
of ~ 10%: The published value from the Hubble Space Telescope Key Project is 71 ± Gkms"-'^ Mpc~^ (Mould, et 
al. 2000). 



III. OBSERVATIONS OF COSMIC MAGNETIC FIELDS 

Observations of galactic and extragalactic magnetic fields can be summarized as follows: 

• Magnetic fields with strength ~ 10/xG are found in spiral galaxies whenever the pertinent observations are made. 
These fields invariably include a large-scale component whose coherence length is comparable to the size of the 
visible disk. There are also small-scale tangled fields with energy densities approximately equal to that of the 
coherent component. 

• The magnetic field of a spiral galaxy often exhibits patterns or symmetries with respect to both the galaxy's 
spin axis and equatorial plane. 

• Magnetic fields are ubiquitous in elliptical galaxies, though in contrast with the fields found in spirals, they 
appear to be random with a coherence length much smaller than the galactic scale. Magnetic fields have also 
been observed in barred and irregular galaxies. 

• Microgauss magnetic fields have been observed in the intracluster medium of a number of rich clusters. The 
coherence length of these fields is comparable to the scale of the cluster galaxies. 

• There is compelling evidence for galactic-scale magnetic fields in a redshift z ~ 0.4 spiral. In addition, microgauss 
fields have been detected in radio galaxies at z > 2. Magnetic fields may also exist in damped Lya systems at 
cosmological redshifts. 
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• There are no detections of purely cosmological fields (i.e., fields not associated with gravitationally bound or 
collapsing structures). Constraints on cosmological magnetic fields have been derived by considering their effect 
on big bang nucleosynthesis, the cosmic microwave background, and polarized radiation from extragalactic 
sources. 

These points will be discussed in detail. Before doing so, we describe the four most common methods used to 
study astrophysical magnetic fields. A more thorough discussion of observational techniques can be found in various 
references including Ruzmaikin, Shukurov, and Sokoloff (1988a). 



A. Observational Methods 

1. Synchrotron Emission 

Synchroton emission, the radiation produced by relativistic electrons spiralling along magnetic field lines, is used to 
study magnetic fields in astrophysical sources ranging from pulsars to supcrclusters. The total synchrotron emission 
from a source provides one of the two primary estimates for the strength of magnetic fields in galaxies and clusters 
while the degree of polarization is an important indicator of the field's uniformity and structure. 

For a single electron in a magnetic field B, the emissivity as a function of frequency v and electron energy E is 

J(y,E) « ' (22) 

where is the component of the magnetic field perpendicular to the line of sight, = vl {E /mc^Y' is the so-called 
critical frequency, vj^ = {eB^_/2'Kmc) is the Larmor frequency, and f{x) is a cut-off function which approaches unity 
for X ^ and vanishes rapidly for a; 3> 1. 

The total synchrotron emission from a given source depends on the energy distribution of electrons, ne{E). A 
commonly used class of models is based on a power-law distribution 



/ E \ — ' 

neiE)dE = jieo i—j dE (23) 

assumed to be valid over some range in energy. The exponent 7 is called the spectral index while the constant 
Ueo = UeiEo) sets the normalization of the distribution. A spectral index 7 ~ 2.6 — 3.0 is typical for spiral galaxies. 
The synchrotron emissivity is ji, = J J{v, E)ne{E)dE. Eq. ( p2[ ) shows that synchrotron emission at frequency v is 

dominated by electrons with energy E ~ rrieC^ (u/i^l)^^^ , i.e., v ~ Vc, so that to a good approximation, we can write 
J{i', E) cx B^i^cS {v — Vc)- For the power-law distribution Eq. ( p3| ) we find 



cx n,oz.(i--)/^Sr-'^/^ . (24) 

Alternatively, we can write the distribution of electrons as a function of Vc'. n{vc) = ne{E)dE / di>c oc j^^/i'cB^. (See 
Leahy 1991 for a more detailed discussion). 

The energy density in relativistic electrons is e^e = / n{E)EdE. Thus, the synchrotron emission spectrum can 
be related to the energy density in relativistic electrons e^e and the strength of the magnetic field (Burbidge 1956; 
Pacholczyk 1970; Leahy 1991). It is standard practice to write the total kinetic energy in particles as efc = (1 -I- k) tre 
where k ~ 100 is a constant (see, for example, Ginzberg & Syrovatskii 1964; Cesarsky 1980). The total energy 
(kinetic plus field) is therefore etot = (1 + fc) Cre + es- One can estimate the magnetic field strength either by 
assuming equipartition ((1 -I- k) tre — ^b) or by minimizing etot with respect to B. 

The standard calculation of e^e uses a fixed integration interval in frequency, vl v < uu: 



ere = / En{vc)dve cx B-^/'^Q^ S^{vo) (25) 

where Su is the total flux density, O is the angular size of the source, and i/q is a characteristic frequency between 
Vl and vjj- Assuming either equipartition or minimum energy, this expression leads to an estimate for B of the form 
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Bcq (X sV Q '^^^ ■ However Beck et al. (1996) and Beck (2000) pointed out that a fixed frequency range corresponds to 
different ranges in energy for different values of the magnetic field (see also Leahy (1991) and references therein). From 
Eq. ( [2^ ) we have n^o oc j{v)v^'^^^^ /'^ B ^^^'^^'^^ ^"^ . Integrating over a fixed energy interval gives e^e Q'^S^B^^^^^^/'^ 
which leads to a minimum energy estimate for B of the form i?cq c>c 5*^/(^+^^9^4/(7+5) 

Interactions between cosmic rays, supernova shock fronts, and magnetic fields can redistribute energy and therefore, 
at some level, the minimum energy condition will be violated. For this reason, the equipartition/minimum energy 
method for estimating the magnetic field strength is under continous debate. Duric (1990) argued that discrepancies of 
more than a factor of 10 between the derived and true values for the magnetic field require rather extreme conditions. 
Essentially, B/B^q sets the scale for the thickness of radio synchrotron halos. A field as small as 0.1i?cq requires higher 
particle energies to explain the synchrotron emission data. However, high energies imply large propagation lengths and 
hence an extended radio halo (scale height ~ 30 kpc) in conflict with observations of typical spiral galaxies. Conversely, 
a field as large as lOBeq would confine particles to a thin disk (~ 300 pc) again in conflict with observations. 

In the Galaxy, the validity of the equipartition assumption can be tested because we have direct measurements of 
the local cosmic-ray electron energy density and independent estimates of the local cosmic-ray proton density from 
diffuse continuum 7-rays. A combination of the radio synchrotron emission measurements with these results yields a 
field strength in excellent agreement with the results of equipartition arguments (Beck 2002). 

While synchroton radiation from a single electron is elliptically polarized, the emission from an ensemble of electrons 
is only partially polarized. The polarization degree p is defined as the ratio of the intensity of linearly polarized 
radiation to the total intensity. For a regular magnetic field and power-law electron distribution (Eq. (p3|)) p is fixed 
by the spectral index 7. In particular, if the source is optically thin with respect to synchrotron emission (a good 
assumption for galaxies and clusters). 



^ = ^^ = 7T773 (2^) 

(Ginzburg & Syrovatskii 1964; Ruzmaikin, Shukurov, & Sokoloff 1988a). For values of 7 appropriate to spiral galaxies, 
this implies a polarization degree in the range p = 0.72 — 0.74. The observed values — p ~ 0.1 — 0.2 for the typical 
spiral — are much smaller. 

There are various effects which can lead to the depolarization of the synchrotron emission observed in spiral 
galaxies. These effects include the presence of a fiuctuating component to the magnetic field, inhomogeneities in the 
magneto-ionic medium and relativistic electron density, Faraday depolarization (see below) and beam-smearing (see, 
for example, Sokoloff et al. 1998). Heuristic arguments by Burn (1966) suggest that for the first of these effects, the 
polarization degree is reduced by a factor equal to the ratio of the energy density of the regular field B to the energy 
density of the total field: 



P = PH^- (27) 

(This expression is useful only in a statistical sense since one does not know a priori the direction of the regular field.) 
Thus, perhaps only ~ 25% of the total magnetic field energy in a typical spiral is associated with the large-scale 
component. Of course, the ratio B/B would be higher if other depolarization effects were important. 



2. Faraday rotation 

Electromagnetic waves, propagating through a region of both magnetic field and free electrons, experience Faraday 
rotation wherein left and right-circular polarization states travel with different phase velocities. For linearly polarized 
radiation, this results in a rotation with time (or equivalently path length) of the electric field vector by an angle 

f = ^^^^^ ne(OS|| {l)dl + ipo (28) 

where rue is the mass of the electron, A is the wavelength of the radiation, ipo is the initial polarization angle, and B\\ 
is the line-of-sight component of the magnetic field. Here, ne{l) is the density of thermal electrons along the line of 
sight from the source (/ — Is) to the observer (/ = 0). tp is usually written in terms of the rotation measure, RM: 
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V = (RM) + ipo (29) 

where 

^ I' -mm 

In general, the polarization angle must be measured at three or more wavelengths in order to determine RM accurately 
and remove the (f = Lp ztmr degeneracy. 

By convention, RM is positive (negative) for a magnetic field directed toward (away from) the observer. The 
Faraday rotation angle includes contributions from all magnetized regions along the line of sight to the source. 
Following Kronberg & Perry (1982) we decompose RM into three basic components: 



RM = RMg + RMs + RMig (31) 

where BMg, RAIg., and RMig are respectively the contributions to the rotation measure due to the Galaxy, the source 
itself, and the intergalactic medium. 

Faraday rotation from an extended source leads to a decrease in the polarization: The combined signal from waves 
originating in different regions of the source will experience different amounts of Faraday rotation thus leading to a 
spread in polarization directions. Faraday depolarization can, in fact, be a useful measure of magnetic field in the 
foreground of a source of polarized synchrotron emission. 



3. Zeeman Splitting 

In vacuum, the electronic energy levels of an atom are independent of the direction of its angular momentum vector. 
A magnetic field lifts this degeneracy by picking out a particular direction in space. If the total angular momentum 
of an atom is J (= spin S plus orbital angular momentum L) there will be 2j + 1 levels where j is the quantum 
number associated with J. The splitting between neighboring levels is AE = gjiB where g is the Lande factor which 
relates the angular momcntmn of an atom to its magnetic moment and /i = eh/2m,,c = 9.3 x 10^^^ erg G""'^ is the 
Bohr magneton. This effect, known as Zeeman splitting, is of historical importance as it was used by Hale (1908) to 
discover magnetic fields in sunspots, providing the first known example of extraterrestrial magnetic fields. 

Zeeman splitting provides the most direct method available for observing astrophysical magnetic fields. Once AE 
is measured, B can be determined without additional assumptions. Moreover, Zeeman splitting is sensitive to the 
regular magnetic field at the source. By contrast, synchrotron emission and Faraday rotation probe the line-of-sight 
magnetic field. 

Unfortunately, the Zeeman effect is extremely difficult to observe. The line shift associated with the energy splitting 

is 



Au , , f B \ /Hz\ 

For the two most common spectral lines in Zeeman-effect observations — the 21 cm line for neutral hydrogen and 
the 18 cm OH line for molecular clouds — Av/v ~ 10~^5f (B/i^G). A shift of this amplitude is to be compared with 

Doppler broadening, Avjv ~ fr/c ~ 6 x 10^^ (T/lOO X)^^^ where vt and T arc the mean thermal velocity and 
temperature of the atoms respectively. Therefore Zeeman splitting is more aptly described as abnormal broadening, 
i.e., a change in shape of a thermally broadened line. Positive detections have been restricted to regions of low 
temperature and high magnetic field. 

Within the Galaxy, Zeeman effect measurements have provided information on the magnetic field in star forming 
regions and near the Galactic center. Of particular interest are studies of Zeeman splitting in water and OH masers. 
Rcid & Silvcrstein (1990), for example, used observations of 17 OB masers to map the large-scale magnetic field 
of the Galaxy. Their results are consistent with those found in radio observations and, as they stress, provide in 
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situ measurements of the magnetic field as opposed to the integrated field along the line-of-sight. Measurements of 
Zecman splitting of the 21 cm line have been carried out for a variety of objects. Kazes, Troland, & Crutcher (1991), 
for example, report positive detections in High Velocity HI clouds as well as the active galaxy NGC 1275 in the 
Perseus cluster. However, Verschuur (1995) has challenged these results, suggesting that the claimed detections are 
spurious signals, the result of confusion between the main beam of the telescope and its sidelobes. Thus, at present, 
there are no confirmed detections of Zeeman splitting in systems beyond the Galaxy. 

4. Polarization of Optical Starlight 

Polarized light from stars can reveal the presence of large-scale magnetic fields in our Galaxy and those nearby. The 
first observations of polarized starlight were made by Hiltner (1949a, b) and Hall (1949). Hiltner was attempting to 
observe polarized radiation produced in the atmosphere of stars by studying eclipsing binary systems. He expected to 
find time-variable polarization levels of 1-2%. Instead, he found polarization levels as high as 10% for some stars but 
not others. While the polarization degree for individual stars did not show the expected time-variability, polarization 
levels appeared to correlate with position in the sky. This observation led to the conjecture that a new property of the 
interstellar medium (ISM) had been discovered. Coincidentally, it was just at this time that Alfvcn (1949) and Fermi 
(1949) were proposing the existence of a galactic magnetic field as a means of confining cosmic rays (See Trimble 1990 
for a further discussion of the early history of this subject). A connection between polarized starlight and a galactic 
magnetic field was made by Davis and Greenstein (1951) who suggested that elongated dust grains would have a 
preferred orientation in a magnetic field: for prolate grains, one of the short axes would coincide with the direction 
of the magnetic field. The grains, in turn, preferentially absorb light polarized along the long axis of the grain, i.e., 
perpendicular to the field. The net result is that the transmitted radiation has a polarization direction parallel to the 
magnetic field. 

Polarization of optical starlight has limited value as a probe of extragalactic magnetic fields for three reasons. First, 
there is at least one other eflFect that can lead to polarization of starlight, namely anisotropic scattering in the ISM. 
Second, the starlight polarization effect is self-obscuring since it depends on extinction. There is approximately one 
magnitude of visual extinction for each 3% of polarization (see, for example, Scarrott, Ward- Thompson, & Warren- 
Smith 1987). In other words, a 10% polarization effect must go hand in hand with a factor of 20 reduction in luminosity. 
Finally, the precise mechanism by which dust grains are oriented in a magnetic field is not well understood (see, for 
example, the review by Lazarian, Goodman, & Myers 1997). 

Polarized starlight does provide information that is complementary to what can be obtained from radio observations. 
The classic polarization study by Mathewson & Ford (1970) of 1800 stars in the Galaxy provides a vivid picture of 
a field that is primarily in the Galactic plane, but with several prominent features rising above and below the plane. 
In addition, there are examples of galaxies where a spiral pattern of polarized optical radiation has been observed 
including NGC 6946 (Fcndt, Beck, & Neininger 1998), M51 (Scarrott, Ward-Thompson, & Warren-Smith 1987), and 
NGC 1068 (Scarrott et al. 1991). The optical polarization map of M51, for example, suggests that its magnetic field 
takes the form of an open spiral which extends from within 200 pc of the galactic center out to at least 5 kpc. Radio 
polarization data also indicates a spiral structure for the magnetic field for this galaxy providing information on 
the magnetic configuration from 3 kpc to 15 kpc (Berkhuijsen et al. 1997). Nevertheless, it is sometimes difficult to 
reconcile the optical and radio data. Over much of the M51 disk, the data indicates that the same magnetic field gives 
rise to radio synchrotron emission and to the alignment of dust grains (Davis-Greenstcin mechanism). However in one 
quadrant of the galaxy, the direction of the derived field lines differ by ^ 60° suggesting that either the magnetic fields 
responsible for the radio and optical polarization reside in different layers of the ISM or that the optical polarization 
is produced by a mechanism other than the alignment of dust grains by the magnetic field. 

B. Spiral Galaxies 

Spiral galaxies are a favorite laboratory for the study of cosmic magnetic fields. There now exist estimates for the 
magnetic field strength in well over 100 spirals and, for a sizable subset of those galaxies, detailed studies of their 
magnetic structure and morphology. 

1. Field Strength 

The magnetic field of the Galaxy has been studied through synchrotron emission, Faraday rotation, optical po- 
larization, and Zeeman splitting. The latter provides a direct determination of the in situ magnetic field at specific 
sites in the Galaxy. Measurements of the 21-cm Zeeman effect in Galactic HI regions reveal regular magnetic fields 
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with i? ~ 2 — 10 yuG, the higher values being found in dark clouds and HI shells (Heiles 1990 and references therein). 
Similar values for the Galactic field have been obtained from Faraday rotation surveys of galactic and cxtragalactic 
sources (i.e., estimates of RMg). Manchester (1974) has compiled RM data for 38 nearby pulsars and was able to 
extract the Galactic contribution. He concluded that the coherent component of the local magnetic field is primarily 
toroidal with a strength B ~ 2.2 ± 0.4 /xG. Subsequent RM studies confirmed this result and provided information 
on the global structure of the Galactic magnetic field (see for example Rand & Lyne 1994 and also Frick, Stepanov, 
Shukurov, & Sokoloff (2001) who describe a new method for analysing RM data based on wavelets). 

Early estimates of the strength of the magnetic field from synchrotron data were derived by Phillipps et al. (1981). 
Their analysis was based on a model for Galactic synchrotron emission in which the magnetic field in the Galaxy is 
decomposed into regular and tangled components. An excellent fit to the data was obtained when each component 
was assumed to have a value of 3 fiG. More recent estimates give 4/iG for the regular and 6/zG for the total local 
field strength (Beck 2002). 

Magnetic fields in other galaxies are studied primarily through synchrotron and Faraday rotation observations. An 
interesting case is provided by M31. Polarized radio emission in this galaxy is confined to a prominent ring ^ lOkpc 
from the galaxy's center. The equipartition field strength in the ring is found to be ^ AjjbG for both regular and 
random components. 

Fitt & Alexander (1993) applied the minimum energy method to a sample of 146 late-type galaxies. The distribution 
of field strengths across the sample was found to be relatively narrow with an average value of (-Beq) — 11 ± 4/xG 
(using k = 100), in agreement with earlier work by Hummel et al. (1988). The magnetic field strength does not appear 
to depend strongly on galaxy type although early-type galaxies have a slightly higher mean. 

A few galaxies have anomalously strong magnetic fields. A favorite example is M82 where the field strength, derived 
from radio continuum observations, is ~ 50 /xG (Klein, Wielebinski, & Morsi 1988). This galaxy is characterized by 
an extraordinarily high star formation rate. 

2. Global Structure of the Magnetic Field in Spirals 

Analysis of RM data as well as polarization maps of synchrotron emission can be used to determine the structure 
of magnetic fields in galaxies. It is common practice to classify the magnetic field configurations in disk galaxies 
according to their symmetry properties under rotations about the spin axis of the galaxy. The simplest examples 
are the axisymmetric and bisymmctric spiral patterns shown in Figure 1. In principle, an RM map can distinguish 
between the different possibilities (Tosa & Fujimoto 1978; Sofue, Fujimoto, & Wielebinski 1986). For example, one can 
plot RM as a function of the azimuthal angle <j) at fixed physical distance from the galactic center. The result will be 
a single (double) periodic distribution for a pure axisymmetric (bisymmetric) field configuration. The RM-0 method 
has a number of weaknesses as outlined in Ruzmaikin, Sokoloff, Shukurov, & Beck (1990) and Sokoloff, Shukurov, 
& Krause (1992). In particular, the method has difficulty disentangling a magnetic field configuration that consists 
of a superposition of different modes. In addition, determination of the RM is plagued by the "riTr degeneracy" and 
therefore observations at a number of wavelengths is required. An alternative is to consider the polarization angle 
'0 as a function of (j). and model as a Fourier series: = a„ cos (n(;i)) + 6„sin(n(/)). The coefficients a„ 

and bn then provide a picture of the azimuthal structure of the field. Of course, if an estimate of the field strength 
is desired, multiwavelength observations are again required (Ruzmaikin, Sokoloff, Shukurov, & Beck 1990; Sokoloff, 
Shukurov, & Krause 1992). 

In M31, both RM{(j)) and "0(0) methods suggest strongly that the regular magnetic field in the outer parts of 
the galaxy (outside the synchroton emission ring) is described well by an axisymmetric field. Inside the ring, the 
field is more complicated and appears to have a significant admixture of either m = 1 or m = 2 modes (Ruzmaikin, 
Sokoloff, Shukurov, & Beck 1990). These higher harmonics may be an indication that the dynamo is modulated by 
the two-arm spiral structure observed in this region of the galaxy. The polarized synchrotron emissivity along the 
ring may provide a further clue as to the structure of the magnetic field. The emissivity is highly asymmetric — in 
general much stronger along the minor axis of the galaxy. Urbanik, Otmianowska-Mazur, & Beck (1994) suggested 
that this pattern in emission are better explained by a superposition of helical flux tubes that wind along the axis 
of the ring rather than a pure azimuthal field. (For a further discussion of helical fiux tubes in the context of the 
aw-dynamo see Donner & Brandenburg 1990). 

Field configurations in disk galaxies can also be classifled according to their symmetry properties with respect to 
reflections about the central plane of the galaxy. Symmetric, or even parity fleld conflgurations are labeled Sm where, 
as before, m is the azimuthal mode number. Antisymmetric or odd parity solutions are labelled Am. Thus as SO field 
configuration is axisymmetric (about the spin axis) and symmetric about the equatorial plane. An AO configuration 
is also axisymmetric but is antisymmetric with respect to the equatorial plane. As shown in Figure 2, the poloidal 
component of a symmetric (antisymmetric) field configuration has a quadruple (dipole) structure. 
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FIG. 1 Axisymmetric and bisymmetric field configurations for disk systems along with the corresponding RM vs. ^ plots. 
Top panels show toroidal field lines near the equatorial plane. Lower panels show RM as a function of azimuthal angle 4> for 
observations at the circle (dotted line) in the corresponding top panel. Note that the pitch angle has the opposite sign for the 
two cases shown. 



The parity of a field configuration in a spiral galaxy is extremely difficult to determine. Indeed, evidence in favor 
of one or the other type of symmetry has been weak at best and generally inconclusive (Krause & Beck 1998). One 
carefully studied galaxy is the Milky Way where the magnetic field has been mapped from the RMs of galactic and 
extragalactic radio sources. An analysis by Han, Manchester, Berkhuijsen, & Beck (1997) of over 500 extragalactic 
objects suggests that the field configuration in the inner regions of the Galaxy is antisymmetric about its midplane. 
On the other hand the analysis by Frick ct al. (2001) indicates that the field in the solar neighborhood is symmetric. 
Evidently, the parity of the field configuration can change from one part of a galaxy to another. A second well-studied 
case is M31 where an analysis of the RM across its disk suggests that the magnetic field is symmetric about the 
equatorial plane, i.e., an even-parity axisymmetric (SO) configuration (Han, Beck, & Berkhuijsen 1998). 

Among SO-type galaxies, there is an additional question as to the direction of the magnetic field, namely whether 
the field is oriented inward toward the center of the galaxy or outward (Krause & Beck 1998). The two possibilities 
can be distinguished by comparing the sign of the RM (as a function of position on the disk) with velocity field data. 
Krause & Beck (1998) point out that in four of five galaxies where the field is believed to be axisymmetric, those fields 
appear to be direc;tcd inward. This result is somewhat surprising given that a magnetic dynamo shows no preference 
for one type of orientation over the other. It would be premature to draw conclusions based on such a small sample. 
Nevertheless, if, as new data becomes available, a preference is found for inward over outward directed fields (or more 
realistically, a preference for galaxies that are in the same region of space to have the same orientation), it would 
reveal a preference in initial conditions and therefore speak directly to the question of seed fields. 



3. Connection with Spiral Structure 

Often, the spiral magnetic structures detected in disk galaxies appear to be closely associated with the material 
spiral arms. A possible connection between magnetic and optical spiral structure was first noticed in observations of 
M83 (Sukumar & Allen 1989), IC 342 and M81 (Krause, Hummel, & Beck 1989a, 1989b). A particularly striking 
example of magnetic spiral structure is found in the galaxy NGC 6946, as shown in Figure 3 (Beck & Hoernes 
1995; Prick et al. 2000). In each case, the map of linearly polarized synchrotron emission shows clear evidence for 
spiral magnetic structures across the galactic disk. The magnetic field in IC 342 appears to be an inwardly-directed 
axisymmetric spiral while the field in M81 is more suggestive of a bisymmetric configuration (Sofue, Takano, & 
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FIG. 2 Field lines for even (top panel) and odd (bottom panel) configurations. Shown are cross-sections perpendicular to the 
equatorial plane and containing the symmetry axis of the galaxy (i.e., poloidal planes). The toroidal field is indicated by an 
'x' (field out of the page) or 'dot' (field into the page). 



Fujimoto 1980; Krause, Hummel, & Beck 1989a, 1989b; Krause 1990). In many cases, magnetic spiral arms are 
strongest in the regions between the optical spiral arms but otherwise share the properties (e.g., pitch angle) of their 
optical counterparts. These observations suggest that either the dynamo is more efficient in the interarm regions or 
that magnetic fields arc disrupted in the material arms. For example, Mcstcl & Subramanian (1991, 1993) proposed 
that the a-efFect of the standard dynamo contains a non-axisymmetric contribution whose configuration is similar to 
that of the material spiral arms. The justification comes from one version of spiral arm theory in which the material 
arm generates a spiral shock in the interstellar gas. The jump in vorticity in the shock may yield an enhanced a- 
effect with a spiral structure. Further theoretical ideas along these lines were developed by Shukurov (1998) and a 
variety of numerical simulations which purport to include nonaxisymmctric turbulence have been able to reproduce 
the magnetic spiral structures found in disk galaxies (Rohde & Elstner 1998; Rohde, Beck, & Elstner 1999; Elstner, 
Otmianowska-Mazur, von Linden, & Urbanik 2000). Along somewhat different lines, Fan & Lou (1996) attempted to 
explain spiral magnetic arms in terms of both slow and fast magnetohydrodynamic waves. 

Recently, Beck et al. (1999) discovered magnetic fields in the barred galaxy NGC 1097. Models of barred galaxies 
predict that gas in the region of the bar is channeled by shocks along highly non-circular orbits. The magnetic field 
in the bar region appears to be aligned with theoretical streamlines suggesting that the field is mostly frozen into the 
gas flow in contrast with what is expected for a dynamo-generated field. The implication is that a dynamo is required 
to generate new field but that inside the bar simple stretching by the gas flow is the dominant process (see Moss et 
al. 2001). 
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FIG. 3 Polarized synchrotron intensity (contours) and magnetic field orientation of NGC 6946 (obtained by rotating iJ-vectors 
by 90°) observed at A6.2cni with the VLA (12.5 arcsec synthesized beam) and combined with extended emission observed 
with the Effelsberg 100 — m telescope (2.5 arcmin resolution). The lengths of the vectors are proportional to the degree of 
polarization. (Prom Beck & Hoernes 1996.) 

4. Halo Fields 

Radio observations of magnetic fields in edge-on spiral galaxies suggest that in most cases the dominant component 
of the magnetic field is parallel to the disk plane (Dumke, Krause, Wielebinski, & Klein 1995). However, for at least 
some galaxies, magnetic fields are found to extend well away from the disk plane and have strong vertical components. 
Hummel, Beck, and Dahlcm (1991) mapped two such galaxies, NGC 4631 and NGC 891, in linearly polarized radio 
emission and found fields with strength ~ 5 and ~ 8 /xG respectively with scale heights ~ 5 — 10 kpc. The fields in these 
two galaxies have rather different characteristics: In NGC 4631 (Figure 4), numerous prominent radio spurs are found 
throughout the halo. In all cases where the magnetic field can be determined, the field follows these spurs (GoUa & 
Hummel 1994). (Recent observations by Tiillmann et al. (2000) revealed similar structures in the edge-on spiral NGC 
5775.) Moreover, the large-scale structure of the field is consistent with that of a dipole configuration (anti-symmetric 
about the galactic plane) as in the bottom panel of Figure 2. The field in NGC 891 is more disorganized, that is, only 
ordered in small regions with no global structure evident. 

Magnetic fields are but one component of the ISM found in the halos of spiral galaxies. Gas (which exists in many 
different phases), stars, cosmic rays, and interstellar dust, are also present. Moreover, the disk and halo couple as 
material flows out from the disk and into the halo only to eventually fall back completing a complex circulation of 
matter (sec, for example, Dahlem 1997). At present, it is not clear whether halo fields arc the result of dynamo 
action in the halo or alternatively, fields produced in the disk and carried into the halo by galactic winds or magnetic 
buoyancy (see Section IV.G). 




FIG. 4 Total radio emission and B vectors of polarized emission of NGC 4631 at A22cm (VLA, 70" synthesized beam). The B 
vectors have been corrected for Faraday rotation; their length is proportional to the polarized intensity (From Krause & Beck, 
unpublished) 



5. Far Infrared-Radio Continuum Correlation 



An observation that may shed hght on the origin and evolution of galactic magnetic fields is the correlation between 
galactic far infrared (FIR) emission and radio continuum emission. This correlation was first discussed by Dickey & 
Salpeter (1984) and de Jong, Klein, Wielebinski, & Wunderlich (1985). It is valid for various types of galaxies including 
spirals, irregulars, and cluster galaxies and has been established for over four orders of magnitude in luminosity (see 
Niklas & Beck (1997) and references therein). The correlation is intriguing because the FIR and radio continuum 
emissions arc so different. The former is thermal and presumably related to the star formation rate (SFR). The latter 
is mostly nonthermal and produced by relativistic electrons in a magnetic field. Various proposals to explain this 
correlation have been proposed (For a review, see Niklas & Beck (1997).) Perhaps the most appealing explanation 
is that both the magnetic field strength and the star formation rate depend strongly on the volume density of cool 
gas (Niklas & Beck 1997). Magnetic field lines arc anchored in gas clouds (Parker 1966) and therefore a high number 
density of clouds implies a high density of magnetic field lines. Likewise, there are strong arguments in favor of a 
correlation between gas density and the SFR of the form SFRoc p" (Schmidt 1959). With an index n = 1.4 ± 0.3, 
taken from survey data of thermal radio emission (assumed to be an indicator of the SFR), where able to provide a 
self-consistent picture of the FIR and radio continuum correlation. 



C. Elliptical and Irregular Galaxies 



Magnetic fields are ubiquitous in elliptical galaxies though they are difficult to observe because of the paucity 
of relativistic electrons. Nevertheless, their presence is revealed through observations of synchrotron emission. In 
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addition, Faraday rotation has been observed in the polarized radio emission of background objects. One example is 
that of a gravitationally lensed quasar where the two quasar images have rotation measures that differ by lOOradm"^ 
(Greenfield. Roberts. & Burke 1985). The conjecture is that light for one of the images passes through a giant cD 
elliptical galaxy whose magnetic field is responsible for the observed Faraday rotation. A more detailed review of 
the observational literature can be found in Moss & Shukurov (1996). These authors stress that while the evidence 
for microgauss fields in ellipticals is strong, there are no positive detections of polarized synchrotron emission or any 
other manifestation of a regular magnetic field. Thus, while the inferred field strengths are comparable to those found 
in spiral galaxies, the coherence scale for these fields is much smaller than the scale of the galaxy itself. 

Recently, magnetic fields were observed in the dwarf irregular galaxy NGC 4449. The mass of this galaxy is an 
order of magnitude lower than that of the typical spiral and shows only weak signs of global rotation. Nevertheless, 
the regular magnetic field is measured to be 6 — 8 G, comparable to that found in spirals (Chyzy et al. 2000). Large 
domains of non-zero Faraday rotation indicate that the regular field is coherent on the scale of the galaxy. This field 
appears to be composed of two distinct components. First, there is a magnetized ring 2.2 kpc in radius in which clear 
evidence for a regular spiral magnetic field is found. This structure is reminiscent of the one found in M31. Second, 
there are radial "fans" - coherent magnetic structures that extend outward from the central star forming region. Both 
of these components may be explained by dynamo action though the latter may also be due to outflows from the 
galactic center which can stretch magnetic field lines. 

D. Galaxy Clusters 

Galaxy clusters are the largest non-linear systems in the Universe. X-ray observations indicate that they are filled 
with a tenuous hot plasma while radio emission and RM data reveal the presence of magnetic fields. Clusters are 
therefore an ideal laboratory to test theories for the origin of extragalactic magnetic fields (see, for example, Kim, 
Tribble, & Kronberg 1991; Tribble 1993). 

Data from the Einstein, ROSAT, Chandra, and XMM-Newton observatories provide a detailed picture of rich 
galaxy clusters. The intracluster medium is filled with a plasma of temperature T ~ 10^ — 10^ K that emits X-rays 
with energies ^ 1 — lOkeV. Rich clusters appear to be in approximate hydrostatic equilibrium with virial velocities 
~ 1000 kms~^ (see, for example, Sarazin 1986). In some cluster cores, the cooling time for the plasma due to the 
observed X-ray emission is short relative to the dynamical time. As the gas cools, it is compressed and flows inward 
under the combined action of gravity and the thermal pressure of the hot outer gas (Fabian, Nulsen, & Canizares 
1984). These cooling flows are found in elliptical galaxies and groups as well as clusters. The primary evidence for 
cooling fiows comes from X-ray observations. In particular, a sharp peak in the X-ray surface brightness distribution 
is taken as evidence for a cooling flow since it implies that the gas density is rising steeply towards the cluster center 
(see, for example, Fabian 1994). 

A small fraction of rich clusters have observable radio halos. Hanisch (1982) examined data from four well- 
documented examples and found that radio-halo clusters share a number of properties — principally, a large ho- 
mogeneous hot intracluster medium and the absence of a central dominant (cD) galaxy. He concluded that radio 
halos are short-lived phenomena, symptoms of a transient state in the lifetime of a cluster. 

Magnetic fields appear to exist in galaxy clusters regardless of whether there is evidence of cooling flows or extended 
radio emission. Taylor, Barton, & Ge (1994), working from the all-sky X-ray sample of galaxies of Edge et al. (1992), 
concluded that over half of all cooling flow clusters have RM > 8GGradm~^ and a significant number have RM > 
2000 radm^^. Furthermore, they found a direct correlation between the cooling flow rate and the observed RM. 
Estimates of the regular magnetic fleld strength for clusters in their sample range from 0.2 — 3 fiG. 

Evidence for magnetic fields in radio-halo clusters is equally strong. Kim et al. (1990) determined the RM for 

18 sources behind the Coma cluster and derived an intracluster field strength of B ^ 2.5 (L/lOkpc) ' /xG where 
L, a model parameter, is the typical scale over which the field reverses direction. Unfortunately, for most clusters, 
there are no more than a few radio sources strong enough to yield RM measurements. To circumvent this problem, 
several authors, beginning with Lawlcr & Dennison (1982), employed a statistical approach by combining data from 
numerous clusters. For example, Kim, Tribble, & Kronberg (1991) used data from ~ 50 clusters (including radio-halo 
and cooling- fiow clusters), to plot the RM of background radio sources as a function of their impact parameter from 
the respective cluster center. The dispersion in RMs rises from a background level of 15 radm~^ to ~ 200 radm~^ 
near the cluster center revealing the presence of magnetic fields in most, if not all, of the clusters in the sample. 
Recently, Clarke, Kronberg, & Bohringer (2001) completed a similar study of 16 "normal" Abell clusters selected to 
be free of widespread cooling flows and strong radio halos. Once again, the dispersion in RM is found to increase 
dramatically at low impact parameters indicating strong (0.1 — 1 fiG) magnetic fields on scales of order 10 kpc. 

Radio emission is of course produced by relativistic electrons spiralling along magnetic field lines. These same 
electrons can Compton scatter CMB photons producing a non-thermal spectrum of X-rays and 7-rays. At high 
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FIG. 5 WSRT map of the Coma cluster of galaxies at 326 MHz from 6 12-h observing sessions. The projected linear of the 
bridge is ~ 1.5h^^ Mpc. The + symbol marks the location of NGC 4839. The peak surface brightness is 9.1 mjy per beam and 
contours are shown at -1,1,2,. ..,9,10,20, ...100,200,400 times 4 mJy per peam. 

energies, these Compton photons can dominate the thermal X-ray emission of the cluster (e.g., Rephaeli 1979). In 
contrast to synchrotron emission, the flux of the Compton X-rays is a decreasing function of B — jc{i^) c< B^^^+'^'^Z^ 
for the power-law electron distribution in Eq. (|2^) — so that an upper limit on the non-thermal X-ray flux translates 
to a lower limit on the magnetic field strength in the cluster. Using this method, Rephaeli & Gruber (1988) found a 
lower limit ^ 10^^ fiG for several Abell clusters in agreement with the positive detections described above. 

E. Extracluster Fields 

There are hints that magnetic fields exist on supercluster scales. Kim et al. (1989) detected faint radio emission in 
the region between the Coma cluster and the cluster Abell 1367. These two clusters are 40 Mpc apart and define the 
plane of the Coma supercluster. Kim et al. (1989) observed a portion of the supercluster plane using the Westerbork 
Synthesis Radio Telescope. They subtracted emission from discrete sources such as extended radio galaxies and found, 
in the residual map, evidence of a 'bridge' in radio emission (Figure 5). The size of the 'bridge' was estimated to be 
1.5/175^ Mpc in projection where /175 is the Hubble constant Hq in units of 75kms~^ Mpc~^. They concluded that 
the 'bridge' was a feature of the magnetic field of the Coma- Abell 1367 supercluster with a strength, derived from 
minimum energy arguments, of 0.2 — 0.6 /iG. 

Indirect evidence of extracluster magnetic fields may exist in radio observations by Ensslin et al. (2001) of the giant 
radio galaxy NGC 315. New images reveal significant asymmetries and peculiarities in this galaxy. These features 
can be attributed to the motion of the galaxy through a cosmological shock wave 10 — 100 times the dimension of a 
typical cluster. Polarization of the radio emission suggests the presence of a very-large scale magnetic field associated 
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FIG. 6 Map of the M81 bisymmetric spiral magnetic field (Krause 1990) projected to have the same linear scale at « = 0.395 
as that of the jet of PKS 1229-021 at the same z. It is shown superposed on the rotation measure distribution of the jet. The 
shaded area shows the region where RM data have been collected and the dashed outline shows the approximate region of the 
total radiation zone. The ridge line of the jet is shown, and the positions of the maxima and minima in RM are shown by 
circled plus and minus signs indicating magnetic field directions toward and away from the observer (from Kronberg, Perry, & 
Zukowski 1992). 

with the shock. 



F. Galactic Magnetic Fields at Intermediate Redshifts 

Evidence of magnetic fields in galaxies at even moderate redshifts poses a serious challenge to the galactic dynamo 
hypothesis since it would imply that there is limited time available for field amplification. At present, the most 
convincing observations of galactic magnetic fields at intermediate redshifts come from RM studies of radio galaxies 
and quasars. Kronberg, Perry, & Zukowski (1992) obtained an RM map of the radio jet associated with the quasar 
PKS 1229-121. This quasar is known to have a prominent absorption feature presumably due to an intervening object 
at z ~ 0.395. (The intervener has not been imaged optically). Observations indicate^ that the RM changes sign along 
the "ridge line" of the jot in a quasi-oscillatory manner. One plausible explanation is that the intervener is a spiral 
galaxy with a bisymmetric magnetic field as illustrated in Figure 6. Alternatively, the field in the intervening galaxy 
might be axisymmetric with reversals along the radial direction. (A configuration of this type has been suggested for 
the Milky Way (see, for example, Pocdz, Shukurov, & Sokoloff 1993).) 

Athreya et al. (1998) studied 15 high redshift (z > 2) radio galaxies at multiple frequencies in polarized radio 
emission and found significant RMs in almost all of them with several of the objects in the sample having RM > 
1000 radm-2. The highest RM in the sample is 6000 radm'^ for the z ~ 2.17 galaxy 1138-262. RMs of this magnitude 
require microgauss fields that are coherent over several kpc. 

G. Cosmological Magnetic Fields 

A truly cosmological magnetic field is one that cannot be associated with collapsing or virialized structures. Cos- 
mological magnetic fields can include those that exist prior to the epoch of galaxy formation as well as those that 
are coherent on scales greater than the scale of the largest known structures in the Universe, i.e., > 50Mpc. In the 
extreme, one can imagine a field that is essentially uniform across our Hubble volume. At present, we do not know 
whether cosmological magnetic fields exist. 

Observations of magnetic fields in the Coma supercluster and in redshift z ~ 2 radio galaxies hint at the existence of 
widespread cosmological fields and lend credence to the hypothesis that primordial fields, amplified by the collapse of 
a protogalaxy (but not necessarily by dynamo action), become the microgauss fields observed in present-day galaxies 
and clusters. (An even bolder proposal is that magnetic fields play an essential role in galaxy formation. See, for 
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example, Wasserman 1978; Kim, Olinto, & Rosner 1996.) The structure associated with the magnetic 'bridge' in the 
Coma supercluster (Kim et al. 1989) is dynamically young so that there has been little time for dynamo processes 
to operate. The observations by Athreya et al. (1998) and, to a lesser extent, Kronberg, Perry, & Zukowski (1992), 
imply that a similar problem exists on galactic scales. Interest in the primordial field hypothesis has also been fueled 
by challenges to the standard dynamo scenario. 

A detection of sufficiently strong cosmological fields would provide tremendous support to the primordial field 
hypothesis and at the same time open a new observational window to the early Universe. Moreover, since very weak 
cosmological fields can act as seeds for the galactic dynamo the discovery of even the tiniest cosmological field would 
help complete the dynamo paradigm. 

For the time being, we must settle for limits on the strength of cosmological fields. Constraints have been derived 
from Faraday rotation studies of high-redshift sources, anisotropy measurements of the CMB, and predictions of light 
element abundances from big bang nucleosynthesis (BBN). 



1. Faraday Rotation due to a Cosmological Field 

Faraday rotation of radio emission from high redshift sources can be used to study cosmological magnetic fields. 
For a source at a cosmological distance Is , the rotation measure is given by the generalization of Eq. ( ^o|) appropriate 
to an expanding Universe: 



rad m ^ Jo \ "^^^ / V / Mpc 

The factor of (1 + z)~^ accounts for the redshift of the electromagnetic waves as they propagate from source to 
observer. We consider the contribution to this integral from cosmological magnetic fields. If the magnetic field and 
electron density are homogeneous across our Hubble volume, an all-sky RM map will have a dipole component (Sofue, 
Fujimoto, Kawabata 1968; Brecher & Blumenthal 1970; Vallee 1975; Kronberg 1977; Kronberg & Simard-Normandin 
1976; Vallee 1990). The amplitude of this effect depends on the evolution of B and ng. The simplest assumption is 
that the comoving magnetic flux and comoving electron number density are constant, i.e., B(z) = iJg (1 + z) and 
ne{z) — rieo (1 + z)"^. The cosmological component of the RM is then 
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where is the angle between the source and the magnetic field. 



F{z) = ^ r dz (1 + zf ^ (35) 
c In dz 



and 



^£ = (1 + Z)-' (an (1 + zf + (1 - an - ^a) (1 + zf + Ua) . (36) 

In Figure 7, we plot F and RAIig as a function of Zs for selected cosmological models. The path length to a source and 
hence the cosmological contribution to the RM are increasing functions of Zg, as is evident in Figure 7. In addition, 
for fixed Zg, the RMig is greater in low-flm models than in the Einstein-de Sitter model, a reflection of the fact that 
the path length per unit redshift interval is greater in those models. 

Eqs. (^)-(|3^), together with RM data for high-redshift galaxies and quasars, can be used to constrain the strength 
of Hubble-scale magnetic fields. The difficulty is that the source and Galaxy contributions to the RM are unknown. 
(Indeed, Sofue, Fujimoto, and Kawabata (1968) reported a positive detection of a 10~^ G cosmological field, a result 
which was refuted by subsequent studies.) By and large, the Galactic contribution is RMg < 200 radm^^ and in 
general decreases with increasing angle relative to the Galactic plane. However, Kronberg and Simard-Normandin 
(1976) found that even at high Galactic latitude, some objects have RM > 200 rad m~^. In particular, the high 
Galactic latitude subsample that they considered was evidently composed of two distinct populations, one with 
{RM'^f/'^ ~ 50radm~^ and another with {RM'^f/'^ ~ 200 rad m~^. (A similar decomposition is not possible at 
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FIG. 7 Intergalactic contribution to the rotation measure [RMig) as a function of source redshift Zs. The left-hand vertical axis 
gives the function F{zs) as defined in Eq. (|35). The right-hand axis gives RMig assuming a 1 fiG field, /175 = 1, rieo = 10~^ cm~^, 
and 6 = 0. Curves shown are for the standard cold dark matter (CDM) model {flm = 1), two LCDM models (CDM with a 
cosmological constant) {Qm = 0.15; Qa = 0.85 and and Sim = 0.3; SIa = 0.7) and an open CDM model {Qm = 0.3; Qa = 0). 



low Galactic latitudes where the contribution to the RM from the Galactic magnetic field is stronger. However, the 
RM sky at low Galactic latitudes does reveal a wealth of structure in the Galactic magnetic field (see, for example, 
Duncan, Reich, Reich, & Fiirst (1999) and Gaensler et al. (2001).) This observation suggests that the best opportunity 
to constrain Hubble-scale magnetic fields comes from the high Galactic latitude, low-RM subsample. For example, 
Vallee (1990) tested for an RM dipole in a sample of 309 galaxies and quasars. The galaxies in this sample extended 
to z ~ 3.6 though most of the objects were at z < 2. Vallee derived an upper limit to RMig of about 2 radm~^, 

corresponding to an upper limit of 6 x 10^^^ G (rieo/lO^^ cm~^) ^ on the strength of the uniform component of a 
cosmological magnetic field. 

If either the electron density or magnetic field vary on scales less than the Hubble distance, the pattern of the 
cosmological contribution to the RM across the sky will be more complicated than a simple dipole. Indeed, if 
variations in and B occur on scales much less than c/Hq, a typical photon from an extragalactic source will pass 
through numerous "Faraday screens" . In this case, the average cosmological RM over the sky will be zero. However, 
cr^j^, the variance of RM, will increase with z. Kronberg and Perry (1982) considered a simple model in which clouds 
of uniform electron density and magnetic field are scattered at random throughout the Universe. The rotation measure 
associated with a single cloud at a redshift z is given by 



__ ^ 2.6(l + z) j [j^,^) (37) 

where rc{z) ^ c/H and Ni, = J ne{z)dl ~ 2ne{z)rc{z) are the cloud radius and electron column density respectively. 
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For simplicity, all clouds at a given redshift are assumed to share the same physical characteristics. The contribution 
to cr|j.jyj from clouds between z and z + dz can be written 

c'4Af(^) = RMrXzf^dz (38) 

where 



-^--rlizWiz)- (39) 

and nc{z) is the number density of clouds in the Universe. Kronberg & Perry (1982) estimated cr|j,j,^ under the 
assumption that the comoving number density of clouds, as well as their comoving size, electron density, and magnetic 
flux are constant, i.e., nc(z) = JT-co (1 + z)'^ , r ^{z) = r {1 -\- z)~^ , ne(2;) = rieo (1 + -z)^ , and (z) = Bo (1 + -z)^. 
The expression for (t|j^j then becomes 
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where F' is an integral similar to the one in Eq. (|35|). 

The Kronberg-Perry model was motivated by spectroscopic observations of QSOs which reveal countless hydrogen 
absorption lines spread out in frequency by the expansion of the Universe. This dense series of lines, known as the Lya- 
forest, implies that there are a large number of neutral hydrogen clouds at cosmological distances. Kronberg & Perry 
(1982) selected model parameters motivated by the Lya-cloud observations of Sargent et al. (1980), specifically, Tieo — 
2.5 X 10~^cm~^, rico — 5Mpc~'^, and — 60 kpc. For these parameters, the estimate for aj^J^.J is disappointingly 
small. For example, in a spatially flat, ilm = 0.3 model, (t^j^^(z = 3) ~ 700 {Bq/ fiG)^ rad^m^"'. Detection above the 
"noise" of the galactic contribution requires ctrm ^ 4 radm^^ or equivalently Bq > 0.1 fiG. A magnetic field of this 
strength would have been well above the equipartition strength for the clouds. 

Clearly, the limits derived from RM data depend on the model one assumes for the Lya-clouds. Recently, Blasi, 
Buries, & Olinto (1999) suggested that the conclusions of Kronberg & Perry (1982) were overly pessimistic. Their 
analysis was motivated by a model for the intergalactic medium by Bi & Davidsen (1997) and Coles & Jones (1991) 
in which the Universe is divided into cells of uniform electron density while the magnetic field is parametrized by 
its coherence length and mean field strength. Random lines of sight are simulated for various model universes. The 
results suggest that a detectable variance in RM is possible for magnetic fields as low as Bq — 6 x 10~^ G. The 
enhanced sensitivity relative to the Kronberg & Perry (1982) observations is primarily due to the larger filling factor 
assumed for the clouds: The clouds in Kronberg & Perry (1982) have a filling factor of order 10^'^ while those of 
Blasi, Buries, and Olinto (1999) have a filling factor of ~ 1. 

Further improvements in the use of Faraday rotation to probe cosmological magnetic fields may be achieved by 
looking for correlations in RM (Kolatt, 1998). The correlation function for the RM from sources with an angular 
separation that is small compared to the angular size of the clouds increases as N^, i.e., {RM1RM2) ~ Nf RM'^, 
where Ni is the average number of clouds along the line of sight (Kolatt 1998). By contrast, cr^yi increases linearly 
with Ni^ i.e., (RJVP) ~ NiRM^. Thus, the signal in a correlation map can be enhanced over the ctrm signal by 
an order of magnitude or more. Moreover, in a correlation map, the "noise" from the Galaxy is reduced. Finally, 
the correlation method can provide information about the power spectrum of cosmological fields and the statistical 
properties of the clouds. 

A different approach was taken by Kronberg & Perry (1982) who argued that since the RMig due to Lya-clouds 
is small, a large observed RMig must be due to either gas intrinsic to the QSO or to a few rare gas clouds (e.g., 
a gaseous galactic halo) along the QSO line of sight. The implication is that the distribution of RMig in a sample 
of QSOs will be highly nongaussian (i.e., large for a subset of QSOs but small for many if not most of the others) 
and correlated statistically with redshift and with the presence of damped Lya-systems. Kronberg & Perry (1982), 
Welter, Perry, & Kronberg (1984) and Wolfe, Lanzetta, & Oren (1992) all reported evidence for these trends in RM 
data from QSO surveys. Wolfe, Lanzetta, & Oren (1992), for example, found that in a sample of 116 QSOs, the 5 with 
known damped Lya-systems had large RMig as compared with 35 (i.e., ~ 30%) of those in the rest of the sample. 
However Perry, Watson, & Kronberg (1993) argued that the case for strong magnetic fields in damped Lya-systems 
is unproven, the most serious problems in the Wolfe, Lanzetta, & Oren (1992) analysis arising from the sparsity and 
heterogeneous nature of the data. In particular, since electron densities can vary by at least an order of magnitude, 
a case by case analysis is required. Moreover, a subsequent study by Oren & Wolfe (1995) of an even larger data set 
found no evidence for magnetic fields in damped Lya-systems. 
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2. Evolution of Magnetic Fields in the Early Universe 



Limits on cosmological magnetic fields from CMB observations and BBN constraints are discussed in the two 
subsections tliat follow. These limits are relevant to models in which magnetic fields arise in the very early Universe 
(see Section V.C). In this subsection, we discuss briefly the pre-recombination evolution of magnetic fields. 

During most of the radiation-dominated era, magnetic fields are frozen into the cosmic plasma. So long as this is 
the case, a magnetic field, coherent on a scale L at a time ti, will evolve, by a later time t2 according to the relation 



Jedamzik, Katalinic, & Olinto (1998) pointed out that at certain epochs in the early Universe, magnetic field energy 
is converted into heat in a process analogous to Silk damping (Silk 1968). In particular, at recombination, the photon 
mean free path and hence radiation diffusion length scale becomes large and magnetic field energy is dissipated. The 
damping of different MHD modes — Alfven, fast magnetosonic, and slow magnetosonic — is a complex problem and 
the interested reader is referred to Jedamzik, Katalinic, & Olinto (1998). In short, modes whose wavelength is larger 
than the Silk damping scale at decoupling (comoving length XgHk — 50 Mpc) are unaffected while damping below the 
Silk scale depends on the type of mode and the strength of the magnetic field. 



3. Limits from CMB Anisotropy Measurements 

A magnetic field, present at decoupling (z^ ~ 1100) and homogeneous on scales larger than the horizon at that 
time, causes the Universe to expand at different rates in different directions. Since anisotropic expansion of this type 
distorts the CMB, measurements of the CMB angular power spectrum imply limits on the cosmological magnetic 
fields (Zel'dovich & Novikov 1983; Madsen 1989; Barrow, Ferreira, & Silk 1997). 

The influence of large-scale magnetic fields on the CMB is easy to understand (see, for example, Madsen 1989). 
Consider a universe that is homogeneous and anisotropic where the isotropy is broken by a magnetic field that is 
unidirectional but spatially homogeneous. Expansion of the spacetime along the direction of the field stretches the 
field lines and must therefore do work against magnetic tension. Conversely, expansion orthogonal to the direction of 
the field is aided by magnetic pressure. Thus, the Universe expands more slowly along the direction of the field and 
hence the cosmological redshift of an object in this direction is reduced relative to what it would be in a universe in 
which B = 0. 

Zel'dovich & Novikov (1983) and Madsen (1989) considered a spatially flat model universe that contains a homo- 
geneous magnetic field. The spacetime of this model is Bianchi type I, the simplest of the nine homogeneous and 
anisotropic three-dimensional metrics known collectively as the Bianchi spacetimes. The analysis is easily extended 
to open homogeneous anisotropic spacetimes (i.e., universes with negative spatial curvature), known as Bianchi type 
V (Barrow, Ferreira, & Silk 1997). In the spatially flat case, the model is described by two dimensionless functions of 
time. The first is the ratio of the energy density in the field relative to the energy density in matter: 



o2 

Q ^ 5 ■ (42) 

The second function is the difference of expansion rates orthogonal to and along the field divided by the Hubble 
parameter. The fact that the angular anisotropy of the CMB is small on all angular scales implies that the both of 
these functions are small. Madsen (1989) found that to a good approximation 

AT 

— > 4Qd (43) 

where the subscript 'd' refers to the decoupling epoch. If we assume that the field is frozen into the plasma, then 
B cx (1 -I- z)^ and Q cx (1 -f z). By definition e„i/ri = SH'^/SnG. The constraint implied by Eq. ( p3[ ) can therefore be 
written 

< 3 X lO^^G . (44) 
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While this expression was derived assuming a pure Bianchi type I model (i.e., homogeneous on scales larger than the 
present day horizon) the main contribution to the limit comes from the expansion rate at decoupling and therefore 
the result should be valid for scales as small as the scale of the horizon at decoupling. 

Barrow, Ferreira, & Silk (1997) carried out a more sophisticated statistical analysis based on the 4-year Cosmic 
Background Explorer (COBE) data for angular anisotropy and derived the following limit for primordial fields that 
are coherent on scales larger than the present horizon: 



Measurements of the CMB angular anisotropy spectrum now extend to scales ^ 500 times smaller than the present- 
day horizon (Balbi et al. 2000; Melchiorri et al. 2000; Pryke et al. 2001). These measurements imply that for fields 
with a comoving coherence length > 10 Mpc their strength, when scaled via Eq. ( ^l| ) to the present epoch, must be 
< 10^^ G (Durrer, Kahniashvili, & Yates 1998; Subramanian & Barrow 1998). Future observations should be able to 
detect or limit magnetic fields on even smaller scales. However, the interpretation of any limit placed on magnetic 
fields below the Silk scale is complicated by the fact that such fields are damped by photon diffusion (Jedamzik, 
Katalinic, & Olinto 1998; Subramanian & Barrow 1998). If a particular MHD mode is efficiently damped prior to 
decoupling, then any limit on its amplitude is essentially useless. 

Jedamzik, Katalinic, & Olinto (2000) pointed out that as MHD modes are damped, they heat the baryon-photon 
fluid. Since this process occurs close to the decoupling epoch, it leads to a distortion of the CMB spectrum. Using data 
from the COBE/FIRAS experiment (Fixen 1996), they derived a limit on the magnetic field strength of _B < 3 x 10^* G 
(scaled to the present epoch) between comoving scales ~ 400 pc and 0.6 Mpc. 

The existence of a magnetic field at decoupling may induce a measurable Faraday rotation in the polarization signal 
of the CMB. Kosowsky & Loeb (1996) showed that a primordial field with strength corresponding to a present-day 
value of 10^^ G induces a 1° rotation at 30 GHz and a strategy to measure this effect in future CMB experiments was 
suggested. 

4. Constraints from Big Bang Nucleosynthesis 

Big bang nucleosynthesis (BBN) provides the earliest quantitative test of the standard cosmological model (see, for 
example, Schramm & Turner 1998; Olive, Steigman, & Walker 2000). BBN took place between 10"^ and Is after 
the Big Bang and is responsible for most of the "'He, '^He, D, and ^Li in the Universe. Numerical calculations yield 
detailed predictions of the abundances of these elements which can be compared to observational data. Over the 
years, discrepancies between theory and observation have come and gone. Nevertheless, at present, BBN must be 
counted as an unqualified success of the Big Bang paradigm. 

Magnetic fields can alter the predictions of BBN. Thus, the success of BBN — specifically, the agreement between 
theoretical predictions and observations of the light element abundances — imply limits on the strength of primordial 
fields. Limits of this type were first proposed by Greenstein (1969) and O'Connell & Matese (1970). They identified the 
two primary effects of magnetic fields on BBN: (i) nuclear reaction rates change in the presence of strong magnetic 
fields, and (ii) the magnetic energy density leads to an increase cosmological expansion rate. During the 1990's, 
detailed calculations were carried out by numerous groups including Cheng, Schramm, & Truran (1994), Grasso & 
Rubenstein (1996), Kernan, Starkman, & Vachaspati (1996), and Cheng, Olinto, Schramm, & Truran (1996). Though 
the results from these groups did not always agree (see, for example, Kernen, Starkman, & Vachaspati (1997). The 
general concensus is that the dominant effect comes from the change in the expansion rate due to magnetic field 
energy. 

The effects of a magnetic field on BBN can be understood in terms of the change they induce in the neutron 
fraction. For example, a magnetic field affects the neutron fraction by altering the electron density of states. In a 
uniform magnetic field, the motion of an electron can be decomposed into linear motion along the direction of the 
field and circular motion in the plane perpendicular to the field. According to the principles of quantum mechanics, 
the energy associated with the circular motion is quantized and the total energy of the particle can be written 



where ris — 1,2,... is the quantum number for the different energy eigenstates known as Landau levels (Landau & 
Lifshitz ) which are important for field strengths B > Be = m1c^/eh = 4.4 x lO'^'^G. The (partial) quantization of 
the electron energy implied by Eq. (^6|) therefore changes the density of states of the electrons which in turn affects 
processes such as neutron decay where it leads to an increase in the decay rate. 



Bcosmic < 5 X 10-9/l75f^'/'G . 



(45) 




(46) 
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If the only effect of the magnetic field was to increase nuclear reaction rates, it would lead to a decrease in the 
number of neutrons at the time of BBN and hence a decrease in the ^iJe abundance. However, a magnetic field 
also contributes to the energy density of the Universe and therefore alters the time-temperature relationship. If the 
correlation length of the field is greater than the horizon scale, the field causes the Universe to expand anisotropically. 
On the other hand, a field whose correlation scale is much smaller than the horizon can be treated as a homogeneous 
and isotropic component of the total energy density of the Universe. In either case, the magnetic field increases the 
overall expansion rate thus decreasing the time over which nucleosynthesis can occur and in particular the time over 
which neutrons can decay. The net result is an increase in the "^He abundance. The helium abundance is fixed when 
the age of the Universe is t ~ 1 s and the temperature is kT ~ 1 MeV. At this time, the energy density of the Universe 
is 2 X lO^^ergcm"^ which is comparable to the energy density in a 6 x 10^^ G magnetic field. The magnetic field 
must be somewhat less than this value so as not to spoil the predictions of BBN. If one assumes that the magnetic 
field scales according to Eq. (^) then this leads to the following constraint on the magnetic field at the present epoch 
B < 10-6 G. 



5. Intergalactic Magnetic Fields and High Energy Cosmic Rays 

Cosmic rays are relativistic particles (primarily electrons and protons with a small admixture of light nuclei and 
antiprotons) that propagate through the Galaxy with energies ranging from 10^ — 10^" eV (Hillas 1998). Their energy 
spectrum is characterized by a power law up to the 'knee' {E ~ 10^^ eV), a slightly steeper power law between the 
knee and the 'ankle', {E ~ 10^^ eV), and a flattened distribution (the ultrahigh energy cosmic rays or UHECRs) 
above the ankle. The origin of the UHECRs is a mystery. Circumstantial evidence suggests that these particles are 
created outside the Galaxy but within 50 — 100 Mpc. Since the gyrosynchrotron radius for a particle in the Galactic 
magnetic field with E > 10^^ eV is larger than the Galaxy, if UHECRs originated in the Galaxy, the arrival direction 
would point back to the source. However, to date, no sources have been identified. On the other hand, protons with 
energies above 5 x 10^^ eV interact with CMB photons producing pions over a mean free path of order 50 — 100 Mpc. 
Therefore if most UHECRs originated at cosmological distances, their energy spectrum would show a distinct drop 
known as the GZK cut-off (Greisen, 1966; Zatsepin & V. A. Kuzmin 1966). The absence of such a cut-off implies that 
UHECRs are produced within 100 Mpc. 

A number of authors have considered the fate of UHECRs that are produced in the local supercluster (LSC) under 
the assumption that the LSC is magnetized (Lemoine, et al. 1997; Blasi & Olinto 1999 and references therein). This 
assumption is reasonable given the detection of magnetic fields in the Coma supercluster (Kim et al. 1989). Blasi & 
Olinto (1999) found that for a LSC field of 10^^ G cosmic rays with energies below 10^^ eV execute a random walk 
as they travel from source to observer while those above 10^° eV follow a relatively straight path. They argued that 
the break in the energy spectrum at the ankle is a consequence of the transition from random walk to free-stream 
propagation. A corollary of this result is that the detection of source counterparts for the 10^^ eV particles (or 
alternatively clustering in the source distribution) would imply a limit on the strength of the magnetic field in the 
LSC. 

The effects of large-scale magnetic fields on UHECRs were also considered by Waxman & Miralda-Escude (1996). 
Following Waxman (1995), Vietri (1995) and Milgrom & Usov (1995) they assumed that the same astrophysical objects 
responsible for 7-ray bursts also produce UHECRs. The expected rate for 7-ray bursts within 100 Mpc is only 1 per 
50 yr. Since Takeda et al. (1998) observed 7 events above 10^" eV over 8 years a dispersion in arrival times of > 50yrs 
for cosmic rays produced in a single burst needs to be invoked. Waxman & Miralda-Escude (1996) proposed that such 
a dispersion is due to the deflections of UHECRs by a large-scale magnetic field. The induced time delay is estimated 
to be T ~ 50 yr (i/10 Mpc) (i?/5 x 10"^°) for a cosmic ray energy of 10^° eV and source at 100 Mpc. Magnetic fields 
have a predictable effect on the angular position and time of flight for cosmic rays of a given energy. Hence, future 
cosmic ray experiments should be able to determine not only whether or not UHECRs are indeed produced by 7-ray 
bursters but whether they are deflected by a large-scale magnetic field on route from source to observer. 

Along rather different lines, Plaga (1995) proposed a technique for detecting cosmological magnetic fields at ex- 
tremely low levels. The idea is to look at the arrival times of 7-ray photons from cosmological sources (e.g., 7-ray 
bursts, flare events in AGN). High energy photons suffer collisions in diffuse extragalactic radiation fields. At 7-ray 
energies, the dominant process is electron-positron pair production. The electrons and positrons can then inverse 
Compton scatter off CMB photons producing high energy photons. An intergalactic magnetic field will deflect the 
electrons and protons and therefore delay the secondary pulse. In principle, this technique could be able to detect 
fields as weak as 10"^^ G! 
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IV. GALACTIC AND EXTRAGALACTIC DYNAMOS 

A magnetic dynamo consists of electrically conducting matter moving in a magnetic field in such a way that the 
induced currents amplify and maintain the original field. The dynamo principle was known in the 1800s though it 
was Larmor (f 9f 9) who first suggested that dynamo processes might be responsible for astrophysical magnetic fields 
such as those found in the Sun and Earth. Steenbeck, Krause, & Radler (1966) recognized the importance of helical 
turbulence for dynamos in stars and planets. Their ideas were soon applied to the problem of galactic magnetic fields 
by Parker (1971) and Vainshtein & Ruzmaikin (1971, 1972). 

Over the years, a standard galactic dynamo model known as the aw-dynamo has emerged whose essential features 
are as follows: Turbulent motions in the ISM driven, for example, by stellar winds, supernova explosions, and 
hydromagnetic instabilities, carry loops of toroidal magnetic field out of the plane of the disk. These loops are twisted 
into the poloidal plane by the Coriolis effect while toriodal field is regenerated from the poloidal field by differential 
rotation. The aw-dynamo can operate in any differentially rotating, turbulent medium and is widely accepted as the 
primary mechanism for the maintenance of magnetic fields in the Sun (Krause & Radler 1980; Zel'dovich, Ruzmaikin, 
& Sokoloff 1983; Ruzmaikin, Sokoloff, & Shukurov 1988a). Its applicability to galaxies has been more controversial 
and numerous variants and alternative models have been proposed. Nevertheless, the general idea that galactic fields 
are maintained by differential rotation and small-scale velocity fluctuations is compelling. More speculative is the 
conjecture that magnetic dynamos act on supergalactic scales. On the other hand, if structure formation proceeds 
hierarchically, it is plausible that dynamo processes operate sequentially from subgalactic to galactic scales. 

An essential feature of a dynamo is its ability to continuously regenerate large-scale magnetic fields. For galaxies, the 
alternative is that magnetic fields are relics of the early Universe. A magnetic field that permeates the protogalactic 
medium will be amplified by compression as a galaxy forms and by differential rotation once the disk is fully developed 
(Hoyle 1958; Piddington 1964, 1972; Kulsrud 1990). The relic field hypothesis has been challenged vigorously by Parker 
(1973b) and others on the grounds that turbulent diffusion destroys a primordial field on a relatively short timescale. 
Since these arguments provide a good introduction to the aw-dynamo, we repeat them below. NB: The dynamo 
hypothesis does not preclude primordial magnetic fields. To the contrary, while a dynamo can amplify existing fields, 
the first fields might have been primordial, i.e., created in the very early Universe. 

After a brief discussion of the primordial field hypothesis, we review the essentials of mean- field dynamo theory 
and the standard aw-dynamo. Particular attention is paid to the assumptions necessary for the development of this 
model. Some of these assumptions have been challenged as being demonstrably false while others are seen as simply 
too restrictive. These concerns have led to alternative models for galactic magnetic fields which are discussed at the 
end of the Section. 



A. Primordial Field Hypothesis 

Soon after the discovery of galactic magnetic fields, Hoyle (1958) began to contemplate their origin. An astrophysical 
battery seemed an implausible explanation for galactic fields since the voltage required to drive the requisite currents, 
V ^ 3x 10^"^ Volts, is so enormous. Instead, Hoyle considered a scenario in which magnetic fields are present ab initio 
in the material that collapses to form a galaxy. Piddington (1964, 1972) championed the primordial field hypothesis 
and developed models for the structure and evolution of an initially homogeneous field (presumed to be of primordial 
origin) in a rotating disk galaxy. The primordial field hypothesis hypothesis has been studied recently by Howard & 
Kulsrud (1997). 

The following idealized example illustrates the difficulties with the primordial field hypothesis (Parker 1973b). 
Consider a differentially rotating disk with angular velocity u = uj{R). (We are using cylindrical (i?, 0, z) coordinates.) 
Suppose that at t = 0, the magnetic field is uniform and lies in the disk plane. Without loss of generality, we can 
orient the x-axis to be along the initial direction of the field, i.e., B(x, 0) = Bqx.. So long as magnetic diffusion and 
backreaction effects are negligible, the field at time t > will be given by 

B {R, cj), t) =B^{h + (cos {uj + cf>t)) (47) 

(see Eq. (|l3|)) where b = b(i?) = cos ujtx + sin ujty. Field lines for an illustrative example are shown in Figure 8. Due 
to differential rotation, the azimuthal component of the field grows linearly with t while gradients in the field grow 
as t^. Note that both the initial and final field configurations are bisymmetric. (Alternatively, if the magnetic field is 
initially oriented along the spin axis of the disk, the field will be axisymmetric but with odd parity about the equatorial 
plane of the disk, that is, an AO configuration (see, for example, Ruzmaikin, Shukurov, & Sokoloff 1988a.) Differential 
rotation, with an axisymmetric velocity field (i.e., lu independent of 0) does not alter the symmetry properties of the 
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FIG. 8 Distortion of magnetic field lines under the action of differential rotation. The upper panel shows the initial homogeneous 
magnetic field configuration. The different line types are for visualization purposes. The lower panel shows the same field lines 
after they have been distorted by differential rotation. 

field. Roughly speaking, we have VB ^ {ujt)^ / L vi^here L is the disk scale length and duj/dr ^ ui/L. Eventually, 
magnetic diffusion becomes important with the diffusion timescale, = B/rj'V^B ^ L"^ / rfixP'f^ , decreasing with t. 

1/3 

Linear growth lasts until t ~ or equivalently t ~ i^L'^ /w'^rj) ' after which the field decays rapidly. In galactic disks, 
L ~ 3kpc, u) ~ 10^^'"' s^^, and rj ~ lO^^cm^s^^ implying a decay time of i ~ 3 x lO^yrs which is much shorter than 
the age of a galaxy. (NB: The timescale for the field to decay depends on the structure of the field. If, for example, 
the field is concentrated in intermittent rope-like structures, the decay time will be much longer than indicated by 
the estimate above (Subramanian 1998).) Equally problematic is the observation that galactic magnetic fields form, 
by and large, a loosely wound spiral as in Figure 3 rather than the tightly wound spiral suggested by Figure 8. The 
implication is that galactic magnetic fields are generated continuously. This argument is similar to the one given 
in support of the hypothesis that spiral arms are (continuously-generated) density waves propagating through the 
galactic disk. 

These conclusions were challenged recently by Howard & Kulsrud (1997) who investigated the primordial field 

hypothesis within the context of a simple model for galaxy formation. In particular, they considered a rotating 
spherical protogalaxy threaded by a constant magnetic field which then collapses to form a disk galaxy. This process 
leads to amplification of the initial field by several orders of magnitude (see also Lesch & Chiba 1995). The field is 
then wound up by differential rotation as described above. Howard & Kulsrud (1997) point out that the magnetic 
field observed in spiral galaxies is an average of the true detailed field. Their contention is that the fine scale structure 
of the field is that of Figure 8 and only appears as a toroidal azimuthal field because of inadequate resolution. 



27 



B. Mean-Field Dynamo Theory 

Most discussions of astrophysical dynamos make use of a mean-field approximation to describe the effects of tur- 
bulence. In addition, backreaction of the field on the fluid is typically ignored so that the evolution of the field 
reduces to a purely kinematic problem. Detailed treatments of mean-field dynamo theory can be found in numerous 
references including Steenbeck, Krause, & Radler (1966), Moffatt (1978), Parker (1979), Krause & Radler (1980), 
Zel'dovich, Ruzmaikin, & Sokoloff 1983, Ruzmaikin, Shukurov, & Sokoloff (1988a), Krause & Wielebinski (1991), 
Beck et al. (1996) and Kulsrud (1999). 

The evolution of a magnetic field in the MHD limit is given by Eq.(P). In a mean-field analysis, we write 



B = B + b V = V + v (48) 

where B and V represent ensemble averages of the magnetic and velocity fields and b and v are the corresponding 
small-scale tangled components. The ensemble average of Eq. (|^) is 



dB 



V X (V X B) + V X (v X b) 



(49) 



with the residual equation 



— = Vx(vxB + Vxb-|-vxb-vxb) . (50) 
ot 

(We have assumed that molecular diffusion can be neglected.) 

The V x b term in Eq. (|50| ) can be eliminated by transforming to the rest frame of the fluid while the v x b terms 
in this equation are usually ignored. Eq. ( |50| ) is then used to eliminate b in favor of v and B in Eq. (^9|), though 
this step requires additional assumptions about the statistical properties of the turbulence. The result, which can be 
derived in a variety of ways (see Ruzmaikin, Shukurov, and Sokoloff 1988a and references therein), is 



a-D 

— = Vx(VxB)+Vxf:. (51) 

£, the effective electromotive force due to turbulent motions of the magnetic fleld as it is carried around by the fluid, 
is often written in terms of two tensors, a and f3: 



£, = a^jBj+(3,,k^ . (52) 

Explicit expressions for the a and /? tensors can be found in various sources including Ruzmaikin, Shukurov, & Sokoloff 
(1988a). 

The classic example of the a-effect is the distortion of a magnetic field line by a localized helical disturbance or 
cyclonic event (Parker 1970). In the context of a galactic dynamo, we can think of a cyclonic event as a plume of 
gas rising above the disk and acted upon by differential rotation and the Coriolis effect. Consider a magnetic field 
configuration that is initially purely toroidal and focus on a single field line at a radius R — Rq. The velocity field of 
the plume is assumed to be constant during the "event" (ti < t <tf) and given by the expression 



Vp(x, t) = — , -, 1 j e V J (53) 

where we have introduced a local Cartesian coordinate system with k = 4> and y = R. By design, V • Vp = 0. The 
field line of interest is initially B(i? — Rq) — Sqx. As before, it is traced by particles that are carried along with the 
fluid. In the absence of rotation, these test particles obey the equation of motion du/dt = (vp • V) Vp where u = u(t) 
is the velocity of a test particle as distinct from the fluid velocity field of the plume Vp. The field line after a cyclonic 
event is shown in Figure 9(a). In a rotating system, it is easiest to follow the evolution of a field line (i.e., the motion 
of the tracer particles) in a frame rotating with angular velocity ijJq = uj(Rq). The equation of motion for the particles 
is then 
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FIG. 9 Cyclonic event as an illustration of the a-effect. We consider a single magnetic field line init ially oriented along the 
X-axis. In (a) we show field line after it has been distorted by the plume velocity field given in Eq. (pSl). In (b) we include 
the Coriolis effect. The dashed line and arrow represent the electric current associated with the loop of magnetic field in the 
yz-plane. 
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^ = (vp • V) Vp - 2a;o X u + . . . (54) 
where the differential rotation term is not shown exphcitly and is ignored in this iUustrative example. The second 



term on the right-hand side of Eq. ( 54 ) describes the Coriolis effect which twists the magnetic plume into the poloidal 
plane (i.e., the yz-plane) as shown in Figure 9(b). Note that the current associated with this loop is anti-parallel to 
the initial magnetic field line. 

Given a specific model for the small-scale velocity field the a and /? tensors can be calculated from first principles 
(see, for example, Ferriere (1992, 1993, 1998)). However, most calculations make use of ad hoc phenomenological 
forms for these functions. A simplifying though questionable assumption is that helical turbulence in galactic disks is 
isotropic. In this case, the a and /3 tensors take the form 



(v.(Vxv)) , 



(55) 



and 



A,fe=/3e„-fc = (56) 

where Sij and tijk are the unit tensor and three-dimensional permutation tensor respectively and r is the correlation 
time of the turbulence. With these expressions for a and (5 the dynamo equation becomes 



— = V X (V X B) + V X (aB - /^V X B) . (57) 

The cyclonic event described above illustrates the dubious nature of the isotropy assumption: Instabilities tend 
to develop in the direction perpendicular to the disk while vorticity is generated along the spin axis of the galaxy. 
Moreover, the largest turbulent eddies in the Galaxy are of order 100 pc in size which is not much smaller than the 
scale height of the disk. Ferriere (1992, 1993, 1998) has calculated the a and /? tensors for a model in which turbulence 
is driven by correlated supernova explosions. The results serve as an explicit example of a case where turbulence is 
decidedly anisotropic. In this model, the a-tensor takes the form 





^sc 





















The off-diagonal terms describe the convection of magnetic field lines away from the equatorial plane (a "V x (V x B)"- 
term in the induction equation) while the diagonal elements characterize the strength of the a-effect along the three 
coordinate axes. However, while the calculations of Ferriere account for the stretching of magnetic field lines by 
the expanding supernova shells and concomitant twisting by the Coriolis effect they do not take into account for 
the a-effect due to turbulence that is undoubtedly generated by the expanding shells. Thus, her analysis probably 
underestimates the strength of the a-effect associated with supernovae. 

Eq. (|5l|) is the starting point for most discussions of astrophysical dynamos. A self-consistent treatment of the fluid 
requires an equation of motion for V which includes a Lorentz force term describing the backreaction of the field on 
the fluid. In addition, backreaction may affect the fluctuation fields v and b thus modifying a and /3. These effects 
will be discussed below. Here, we assume that backreaction is negligible so that V, a, and /? can be specified as model 
inputs that are independent of B. If we further suppose that they are time- independent, the solutions to Eq. ( [5l| ) 
will be of the form B cx expFi where the eigenvalue F depends on boundary conditions for the field. The solutions of 
interest are, of course, ones where F is real and positive. 



C. Disk Dynamos 

We consider an axisymmetric differentially-rotating disk galaxy with a large-scale velocity field V(i?) = a; x R 
where ui — lu{R)z. In general a and (3 are functions of R and z with the proviso that a{~z) — ^a{z). Rotational 
invariance of the dynamo equations suggests the following Fourier decomposition: 
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B(r, t) = ^B,„(r, z, t)e™^. (59) 

m 

One might guess, and detailed calculations confirm, that the growth rate decreases with increasing m so that the fastest 
growing mode has azimuthal symmetry (m = 0) (Ruzmaikin, Shukurov, and SokolofF 1988a). The m = 1 bisymmetric 
mode can also be important, though in general this mode is difficult to excite without introducing nonaxisymmetric 
forms for a, /3, and/or w. We will return to this point below and proceed with a discussion of axisymmetric solutions. 
In component notation, Eq. (|5^) becomes 

^ = -^("^0)+/5(V^5)h (60) 



where the overbar and subscript m{= 0) are omitted for the sake of clarity. It is generally assumed that the a-terms 
in the equation for are small compared to the w-term (i_^., the toroidal field is generated through the action 
of differential rotation rather than turbulence). Eqs. (|60|)-(p2[) are symmetric under the transformation z — + ~z. 
Therefore, the solutions will have definite parity with both even and odd parity solutions (B''" and B~ respectively) 
possible: 



B%(-z) = ±B^R{z) B^(-z)^±B^cj,{z) Bt{-z) = TBt(z) (63) 

For m — Q, B+ and B" correspond respectively to the even parity (5*0) and odd parity (AO) configurations described 
in Section III.B.2 and shown in Figure 2. 

Since the dynamo equations are linear in the fields they admit solutions of the form B cx exp Tt. Furthermore, for 
disk-like geometries, variations in the field with respect to z will be much greater than those with respect to i? — in 
our own Galaxy, the scale height of the disk is a factor of 10 smaller than the scale radius of the disk. This situation 
suggests a quasi-separation of variables of the form 



B {R, z, t) = Q^{R)B{R, z)e^"* (64) 

(see, for example, Ruzmaikin, Shukurov, and Sokoloff 1988a) where n labels the different eigenfunctions Qn and 
eigenvalues r„. (For a more rigorous discussion of the thin-disk limit in which the Schrodinger-type radial equation 
(see below) appears as an approximation to the full (integro-differential) see Priklonsky et al. 2000.) B depends on 
R parametrically, i.e., no derivatives with respect to R appear in the equation for B. With this ansatz, the dynamo 
equations become 

^(i?)B«^-^(aB,)+/3^, (65) 



and 



duj 
d\nR' 



:Bf 



■f3 



dz'^ 



(66) 



R^dRKR dR ^^"^ 



+ (7(i?) - r„) g„ = . 



(67) 
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(Since does not appear in the equations for Br and B^ we may ignore it at this time. Once Br and B^ have been 
determined, Bz may be found from the condition V • B = or alternatively by solving the B^ equation.) 

Eqs. (^)-(|6^) can be solved once a complete set of boundary conditions are specified. For simplicity, we assume 
that the disk is defined by sharp boundaries at z = zth and R = with force-free fields (i.e., V x B = 0) outside 
the disk. These conditions imply that ^^(z — h) = and B]^{z = /i) ~ 0. {Bji{z = h) is identically zero in the limit 
h/Ro — > but may be small for a finite disk.) 

Eqs. (^) and (|6|) make clear the essential features of the standard acj-dynamo. Turbulence via the a-effcct 
generates B/j (and Bz) from B^ while differential rotation regenerates B^. The sequence of steps for an AO dynamo is 
shown schematically in Figure 10. We begin in Figure 10(a) with a pure poloidal dipole-like field. Differential rotation 
stretches field lines creating toroidal field. This is illustrated in Figure 10(b) where a single field line from Figure 10(a) 
is shown after it has been acted upon by differential rotation. The field in the equatorial plane is characterized by 
strong gradients and high magnetic tension. This tension can be relieved either by turbulent diffusion, via the /3-effect, 
or by some other process (e.g., magnetic reconnection, see below). The net result is to decouple the toroidal field in 
the upper and lower hemispheres, as shown in Figure 10(c). Next we assume that cyclonic events occur throughout 
the disk. The toroidal field is distorted in the vertical direction as in Figure 10(d). The loops of vertical field are 
then twisted into the poloidal plane by the Coriolis effect (Figure 10(e), see, also Figure 9(b)). Once again, some 
form of diffusion or dissipation is needed to eliminate magnetic field near the equatorial plane. Provided this occurs, 
poloidal loops in upper and lower hemispheres can combine to yield a dipole-like field which reinforces the original 
field (Figure 10(f)). 

This example illustrates the importance of diffusion for the dynamo. The a and uj effects twist, shear, and stretch 
magnetic field lines but do not create new ones. While they can increase the magnetic field energy in the system, they 
cannot change the net flux through a surface that encloses it. Diffusion eliminates unwanted flux. In the odd parity 
dynamo of Figure 10, this process occurs in the equatorial region. In an even parity dynamo, diffusion allows flux of 
the wrong sign to escape by moving to high galactic latitudes. 

The importance of diffusion can be illustrated by considering Fr, the radial flux through the surface at R ~ Rd 
and F^, the azimuthal flux through a surface at fixed (see, for example, Zel'dovich, Ruzmaikin, & Sokoloff 1983): 



Fr, 



2ir 



Rudcj) 



dz Bf 



F,, 



Rd 



dR 



dz Ba 



(68) 



From Eq. (60) we find 



— — = 2'kRdP^— 
dt oz 



(69) 



where contributions from the a-term in Eq. (^0|) vanish because a{z = 0) = and ^^(z — h) — 0. With (3 = 0, Fr is 
constant, i.e., growing mode solutions require /? 7^ 0. 

The equation for F^ is also interesting (Zel'dovich, Ruzmaikin, & Sokoloff 1983). From Eq. ( |6l] ) we find 



'JA^r^R rdz^RBn+p'''^ 



dt 



dR 



dz 



(70) 



For the first term, we integrate by parts and use the condition V • B = to find 



J dRj dz—RBR = J dRivBz\^. (71) 

Thus, with P = 0, F^ grows linearly with time, but only if field lines exit the disk through the surface at z = ft,. If 
the field lines are confined to the disk, than F^ remains constant even though the energy in the field is increasing due 
to differential rotation. 

Thus, magnetic fiux must be expelled from a galaxy inorder for the net large-scale field to grow. Flux expulsion can 
occur by a number of mechanisms including magnetic buoyance and supernovae or superbubble explosions. The latter 
was considered by Rafikov & Kulsrud (2000) who concluded that the gravitational field of the disk would severely 
limit flux expulsion. Clearly, this issue deserves further attention, perhaps through three-dimensional numerical 
simulations, which can model a galaxy and its immediate environment. 
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Eqs. ( |65| ) and ( |66| ) can be simplified by introducing the dimensionless quantities i = t/Tn, 7 = ITd, ^ = z/h, 
a = a/ao and uj = — {duj/d\nR) /uj^ where td = h? / (3 is the (turbulent) diffusion timcscale for the large-scale field. 
We have 

—jBR^-R^-i&B^) (72) 

7 - ^) = -RuBr (73) 

where we have introduced the dimensionless dynamo coefhcients Ra — aoTo/h and R^^ ~ ^^oTd- Upon rescaling 
the ratio Br/B^ by Ra, we see that Eqs. ( [72[ ) and ( [73| ) can be characterized by the single dimensionless parameter, 
D = RaRuj = ^oOio/r'^h, known as the dynamo number. 



D. Growth Rate for the Galactic Magnetic Field 

Eqs. dfl) and ^ can be solved using standard techniques. In the thin-disk limit, the fastest growing solution has 
positive parity (i.e., the quadrupole (SO) mode). Ruzmaikin, Shukurov, Sokoloff (1988a) have determined the growth 
rate (i.e., eigenvalue 7) in the thin-disk limit as a function of the dynamo number D for various forms of a{z) (e.g., 
a = sinTTz, z, and 6(2) — 0(— z) where Q is the Heaviside function). Growing mode solutions are found for values 
oi D > Dcr where Dcr, the critical dynamo number, varies from 6 to 11 for the different forms of a mentioned above. 
For D 3> -Dcr, 7 oc D^/^. For example, in the case a — sIuttz, the growth rate can be approximated by the fitting 
formula 7 ~ Q.l^D^/^ - 0.23 (Field 1995). 

The actual growth rate depends on the parameters wq, ao, and Id. It is relatively straightforward to determine 
Wo- For example, at the position of the Sun in the Galaxy, w(i? = Rs) = loq = 2A ~ 29kms~^ kpc~^ where A is the 
first Oort constant (Binney and Merrifield, 1998). ao and to are more difficult to estimate. However, if we write 7 
in units of ujq, 



^ - 0.78 I I - 0.23 ' 



LOO \huJo J \ujqtd 



Ra 



1/2 



< 0.78 ^-fj , (74) 

we see that the dependence on td is relatively weak, and in any case one can ignore the second term and derive an 
upper limit for the growth rate or a lower bound on the e-folding time for the field. Given an observed value for the 
field at a time tf (e.g., in the Galaxy, Bf = 3/xG with tf = to) and an estimate for an initial time ti when the dynamo 
begins to operate, one can derive a lower bound for the requisite seed field. 



E. Criticisms of Mean-Field Dynamo Theory 

The validity of mean-field dynamo theory was questioned almost as soon as it was proposed. The most serious 
criticisms revolve around the assertion that backreaction of the magnetic field on the fluid is unimportant, an assump- 
tion that has been challenged on the following grounds: In a highly conducting turbulent fluid, the magnetic fleld on 
small scales builds up rapidly via small-scale dynamo action and also the tangling of weak large-scale fields if they 
are present. In a mean-field dynamo, amplification of the regular field (length scale much larger than the outer scale 
of the turbulent velocity field) takes place on a timescale much longer than the eddy diffusion time associated with 
the turbulence. The difficulty is that Lorentz forces on small scales can react back on the fiuid altering the turbulent 
motions (Piddington 1964, 1972; Cattaneo & Vainshtein 1991; Vainshtein & Cattaneo 1992; Kulsrud & Anderson 
1992). If turbulent motions are suppressed, then so too will be turbulent diffusion and the a-effect effectively shutting 
off the dynamo. 

The crux of the putative problem lies in the high magnetic Reynolds number (or alternatively, high electric 
conductivity) found in astrophysical fluids. In a highly conducting turbulent fluid, magnetic flux tubes are con- 
tinuously stretched into intricate ribbons thus leading to a build up of energy to small scales until backreaction 
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effects become important. Various authors have suggested that the suppression of the a effect takes the form 

a ~ aT / ^1 + (J5/_Beq)^ where p is a constant of order unity and ax — vtI is the standard result calculated in 

the absence of backreaction, vt and I are the characteristic velocity and length scale associated with the turbulence, 
and i?^q — Anpv'^ is the corresponding equipartition field strength (see, for example, Cattaneo & Vainshtein 1991; 
Vainshtein & Cattaneo 1992; Gruzinov & Diamond (1994)). Thus, if p ~ 1, as has often been suggested, dynamo 
action in galaxies is strongly suppressed. 

The form for a-suppression given above can be traced to Zcl'dovich (1957) who argued that the ratio of the rms field 

strength -Brms (principally tangled fields) to the large-scale field strength B is given hy B/ -Brms — Rm ■ This result 
is valid for a two-dimensional turbulent flow and can be understood as follows (Parker 1979; Field 1995): Consider 
a flux tube with cross-sectional width S so that the timescale for Ohmic diffusion is td ~ S^/v ~ {S'^Rm/vtL). 
Suppose that initially there is a regular field B with coherence length L. In a turbulent flow and in the absence of 
backreaction, the length I of a flux tube grows as I ~ Lexp (i/2r2^) where tt ~ Ljvj- is the characteristic turnover 
time for the turbulence. The field strength grows by the same factor (see Eq. (O) with the assumption that the 
flow is incompressible) so that Bj-^^/B ~ [l/L) — exp(t/2TT). Since the flux through the tube is constant, the 
cross-sectional area must decrease by the corresponding amount. In a two-dimensional turbulent flow, the tube 
stretches into a thin ribbon: One dimension shrinks exponentially while the other dimension remains approximately 
constant. The diffusion time, which is governed by the small dimension of the tube, therefore drops exponentially: 
td — {LRm / Vt) Q'i^'P {—t / tt) — ttRm {B/B,-^s) ■ Again, if backreaction is ignored, the system rapidly approaches 
the (Ohmic) diffusive limit when td — tt or equivalently B/By^w — Rm ■ 

1/2 

Of course, once the field reaches its equipartition value B^q = (47rpw|,) , backreaction effects become important. 
For B > Bcq, Tjj is approximately equal to the standard value for Ohmic diffusion, td ~ L'^/v ~ tt/Rm- If 

1/2 

Beq > B > Rj^j Beq, the cascadc to small scales is blocked by backreaction. In this case the diffusion time (assumed 
to be equivalent to the dissipation time tu) is modified by a factor of B^^/RmB < 1 (Cattaneo & Vainshtein 1991). 

— —1/2 

Conversely, for B < Rj^j Bcq, the kinematic assumption applies and td — tt- These results are summarized by the 
phenomenological formula 



rn=rT(l+ ^ % ]. (75) 

V (i?cq/s) +i; 

Two-dimensional simulations by Cattaneo & Vainshtein (1991) appear to support these conclusions. The implication 
is that the (3 coefficient in the mean-field dynamo equation is reduced relative to the canonical turbulent diffusion 

]^ /2 

value, Pt — vL/3, whenever the large-scale field is greater than i?^^^ ' Bcq. 



/3eff = /3t f 1 + ^^^2 ] ■ (76) 

V (Bcq/B) +lj 

Vainshtein & Cattaneo (1992) argued that a is suppressed in a similar manner. In present-day galaxies, B ~ Bcq and 
Rm — 10^*^ implying that f3 and a are reduced by enormous factors. 

Can the critisms described above be reconciled with the standard dynamo hypothesis? A somewhat extreme 
viewpoint, expressed by Kulsrud & Anderson (1992), is that the physical basis for the mean-field dynamo theory 
is invalid and therefore the galactic magnetic field is of primordial origin. While this represents a departure from 
mainstream thought on the problem of galactic magnetic fields, it deserves careful consideration, especially in light of 
observations of microgauss fields in high redshift galaxies and on supercluster scales in the low redshift Universe. 

A second possibility is that the arguments presented against the mean field dynamo contain fundamental flaws. 
Blackman & Field (1999) pointed out that some of the analytic approaches (e.g., Gruzinov & Diamond 1994) do 
not distinguish between turbulent quantities of the "zeroth-order" state (i.e., the state with no large-scale field) and 
higher order quantities. The former are homogeneous, isotropic, and stationary. The same cannot be said of the 
latter. Blackman & Field (1999) are quick to point out that they do not prove that a mean-field dynamo evades the 
backreaction problem, but only that some of the objections can be challenged. 

In a subsequent paper, Blackman & Field (2000) argued that the a-suppression seen in numerical simulations by 
Cattaneo & Vainshtein (1991), Vainshtein & Cattaneo (1992), and Cattaneo & Hughes (1996) are due to imposed 
(periodic) boundary conditions. Indeed, Blackman & Field (2000) noted that no such suppression is seen in the sim- 
ulations of Brandenburg & Donner (1997) where dynamo action in an accretion disk (and without periodic boundary 
conditions imposed) is considered (see below). 
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Subramanian (1998) proposed a different resolution to the backreaction controversy in suggesting that the fields 
generated by small-scale dynamo action do not fill space but rather are concentrated into intermittent rope-like struc- 
tures. In a weakly ionized gas, the small-scale dynamo saturates when the field in the flux tube reaches equipartition 
strength. However, because of the low filling factor, the average energy density of the field is still well below that of 
the turbulent gas and therefore the gas is unaffected. Fields do tend to straighten out on small scales. Lorentz forces 
act on the charged component of the gas which leads to slippage between the ions and magnetic field on one hand, 
and the neutral gas component on the other. This process is known as ambipolar diffusion and has been considered 
by various authors in the context of galactic dynamo theory (e.g., Zweibel 1988; Brandenburg & Zweibel 1994). 

A final possibility is that formulation of the mean-field dynamo in terms of a turbulent eddy diffusion parameter 
/? and turbulent a-effect is in fact wrong but that some other mechanism provides a and /?-like effects which enter 
the mean field dynamo equation more or less as in Eq. (|5^). This point of view was advocated by Parker (1992) who 
proposed a dynamo based on buoyancy of magnetic flux tubes and neutral point reconncction. Parker's model is 
discussed below. 

One final note: Even if the problems at small scales are resolved, backreaction of the field on the fiuid certainly 
becomes important once the energy density in the total magnetic field becomes comparable to the kinetic energy 
associated with the turbulent eddies. In this case, a quasisteady-state "equipartition" system should emerge. The 
apparent approximate equality of magnetic field energy and "turbulent fluid energy" in the Galaxy suggests that 
mature spiral galaxies have already reached equipartition. Observations of microgauss fields in galaxies at cosmological 
redshifts lend further support to this conjecture. Quenching of the dynamo at equipartition field strengths has been 
discussed by a number of authors including Moffatt (1978) and Krause and Radler (1980). A phenomenological 
approach is to replace a in the dynamo equation with 



aB = aT(l + ^) . (77) 

Quenching, via Eq. (|77|) , can be easily incorporated into numerical simulations of mean-field dynamo action (see 
below) . 



F. Numerical Simulations of Disk Dynamos 



Numerical simulations provide an approach to the study of magnetic dynamos that is complementary to the eigen- 
value analysis described above. Fewer assumptions are required and one can include additional effects not easily 
treated analytically. For example: 

• The separation-of- variables ansatz, Eq. (^), which follows from the thin-disk approximation, is no longer nec- 
essary. Simulations allow one to explore a family of models that interpolate between disklike and spherical 
systems. 

• The full mean-field dynamo equation can be used. In particular, the a-term in the toroidal equation is retained 
allowing for the possibility of a so-called a^-dynamo where toroidal as well as poloidal fields are generated by 
small-scale turbulent motions. 

• Nonlinear feedback of the field on the fiuid can be included, albeit in an ad hoc manner, by requiring that a 
and /3 approach zero as B rises above the equipartition strength. Feedback of this nature can be implemented 
by replacing a with the form given in Eq. ([77[). 

• The assumption that turbulence is isotropic can be dropped. In particular, a and (3 can assume forms more 
general than those given in Eq. ( p5| ) and Eq. ([5^). 

• Additional physics, such as the influence density waves and feedback from star formation, can be included. 

• Simulations can be tailored to individual galaxies and can be used to generate synthetic observations such as 
radio continuum, polarization, and RM maps. 

• Ultimately, we may be able to bypass the mean-field approximation by simulating the evolution of small-scale 
magnetic and velocity fields explicitly. 
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To be sure, many of the effects listed above have been addressed through analysis. For example, Mestel & Subra- 
manian (1991) investigated the influence of spiral density waves via a non-axisymmetric a tensor in an attempt to 
explain the existence of bisymmetric magnetic fields in at least some spiral galaxies. Likewise, nonlinear effects have 
been treated analytically by various authors such as Belyanin, Sokoloff, & Shukurov (1994). 

Over the years, galactic dynamo simulations have improved steadily both in terms of the physics that is included 
and the dynamic range achieved. In early investigations (e.g., Elstner, Meinel, & Riidiger 1990) a, and lu were 
assumed to be axisymmetric. Moreover, quenching of a and (3 due to backreaction of the field on the fluid was ignored. 
With these assumptions, the different terms in the series expansion given by Eq. (^9|) decouple. Their evolution can 
be studied numerically by solving the finite-diff'erence analog of Eqs. (|6^)-(|6^) on a two-dimensional Eulerian mesh. 
Simulations of this type have been carried out by Elstner, Meinel, & Beck (1992) who considered a suite of disklike 
galaxy models characterized by an angular velocity profile uj{R) and analytic forms for the a and /^-effect tensors. In 
almost all cases, the preferred configuration was SO. Only when the a-effect was confined to the inner part of the disk 
was the AO configuration dominant. Furthermore, no examples were found where a bisymmetric mode was preferred, 
though in some cases a bisymmetric mode could be excited. 

Brandenburg et al. (1992) used two-dimensional calculations to explore the difference between spherical and disk 
dynamo models. In particular, they considered models in which a disk is embedded in a gaseous halo. Not surprisingly, 
they found that for certain regions of parameter space, mixed parity modes could be excited. The relevance of these 
results to galactic magnetic fields will be discussed below. 

Only with 3D simulations have the true benefits of numerical simulations been realized. Deviations from azimuthal 
symmetry (for example, due to spiral density waves) can be introduced by allowing cj, a, and /? to depend on (p. 
In 3D simulations, one can study mode coupling and in particular, the effects of quenching via Eq. (|77|). Recently 
simulations have included the effects of gas and stars using N-body and hydrodynamic techniques. These simulations 
are able to treat self-consistently the interplay between spiral density waves and magnetic spiral structures (see, for 
example, Elstner et al. 2000). Present day simulations reveal a wide range of magnetic configurations, a reassuring 
result given the rich variety of magnetic structures observed in actual galaxies (see, for example, Panesar & Nelson 
1992; Rohde & Elstner 1998). 

In the simulations discussed so far, the effects of turbulence are treated via the mean field quantities a and /?. High 
performance computers make it possible to study dynamo processes from first principles. The challenge, of course, is 
to achieve the dynamic range necessary to follow small-scale fluctuations in the velocity and magnetic fields. Here, 
we mention several techniques used to simulate MHD phenomena that may soon yield useful results in the study of 
dynamo theory. Eulerian finite-difference codes such as ZEUS (Stone & Norman 1992a, 1992b) are designed to solve 
the equations of ideal MHD. A key feature of these codes is that V • B = is guaranteed explicitly at each timestep. 
Recently, Roettighr, Stone, & Burns (1999) used ZEUS to study the evolution of magnetic fields in merging galaxy 
clusters. They found that the field strength and structure are dramatically altered during the merger. In the early 
stages of the merger, the field is stretched by bulk flows and compressed by shocks. However, significant amplification 
of the field does not occur until the later stages of the merger when gas motions become turbulent. 

Korpi et al. (1999) developed a numerical model for the ISM that includes the effects of rotation, supernova heating, 
and nonideal MHD phenomena. In principle, simulations carried out within the context of this model should be able 
to follow magnetic fleld amplification by dynamo action. However, the calculations performed thusfar have not been 
run over an adequate time span to determine whether or not magnetic field generation by dynamo action has occurred 
(see, also, Shukurov (1999)). 

Today, most cosmological simulations use a Lagrangian or particle description for matter. Stars and dark matter 
are modeled as collisionless particles while gas is treated using smooth particle hydrodynamics (SPH; Monaghan 1992 
and references therein). In SPH, any local physical quantity /(r) (e.g., density, pressure) can be approximated by 
taking a weighted average of the form 




(78) 



where rrii is the mass of the i'th particle. 




(79) 



is the density, and W (r — r^, h) is a user-supplied window function whose characteristic width is h. 

The SPH prescription can be used to model magnetic fields (see, for example, Monaghan 1992). Each "gas particle" 
carries with it an additional physical quantity, a B-vector, whose evolution is followed by solving the SPH analog of 
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the induction equation. In general, particle simulations such as those that use SPH have more dynamic range than 
simulations done on a mesh. One potentially serious problem with using SPH to study MHD phenomena is that 
V • B = is not constrained to vanish. 

Dolag, Bartelmann, & Lesch (1999) used SPH to study the evolution of magnetic fields in a cosmological setting. 
In particular, they followed the formation of a galaxy cluster from the linear regime (rcdshift Zi ~ 15) to the present 
epoch assuming an initial magnetic field Bi ~ 10^^ G. The magnetic field is amplified by a factor ^ 1000, roughly an 
order of magnitude more than what would be expected from simple collapse calculations. The simulations were used 
to generate synthetic RM maps for the simulated clusters which are in good agreement with those observed by Kim 
et al. (1990) and Kim, Kronberg, & Tribble (1991). 

Mean-field dynamo theory is based on the assumption that the coherent component of a magnetic field fills space 
and described properly by field equations. (This assumption is the basis for both the Eulerian and SPH methods 
described above.) An alternative picture, supported to some extent by both simulations and observations, is that the 
magnetic field is confined to flux tubes. For example, Roettighr, Stone, & Burns (1999) found filamentary structures 
in their cluster simulations. Moreover, it may be that magnetic fields in extragalactic objects are filamentary but 
appear smooth because of lack of resolution (Hanasz & Lesch 1993). In fact, filamentary structures have been observed 
in the halo of M82 (Renter et al. 1992). The radio spurs seen in NGC 4631 (GoUa & Hummel 1994) and NGC 5775 
(Tiillman (2000) may also be indicative of individual flux tubes. In addition, the magnetic field in the ISM appears 
to be highly nonuniform, again suggesting that the magnetic field in the ISM may be "ropy" (Heiles 1987). 

The flux tube picture suggests an alternative treatment of astrophysical magnetic fields based on the thin-tube 
approximation (Spruit 1981). The working assumption is that the flux tube radius is much smaller than the charac- 
teristic scales over which physical quantities (e.g., gravitational field, gas density) vary. In this limit all variables, such 
as the field strength and gas density inside the tube, can be treated as functions of their position along the tube. These 
quantities evolve according to a set of (Lagrangian) equations of motion. Eqs. ( [lO| ) and ( [ll| ) are examples of this, 
applicable in the kinematic regime. Thus, the illustrative examples in Figures 8-10 may be regarded as (kinematic) 
flux tube simulations. A recent example of a flux-tube simulation can be found in Vainshtein et al. (1996) where the 
so-called stretch-twist-fold dynamo (Vainshtein & Zel'dovich 1972) was investigated using precisely these equations. 

While the flux-tube model has been used extensively to study solar magnetic fields, its application to the problem 
of galactic and extragalactic magnetic fields has been rather limited. One example is the model of Hanasz & Lesch 
(1993, 1997, 1998) who considered it in the context of the Parker instability (see below). 

G. Diversity in Galactic IVIagnetic Fields 

The variety of field configurations observed in spiral galaxies presents a challenge to the dynamo hypothesis. In 
a thin-disk aw-dynamo, the fastest growing mode has SO symmetry and therefore, except in cases of very special 
initial conditions, this mode dominates at late times. However, there are a number of edge-on spiral galaxies in which 
significant poloidal fields have been observed. For example, polarization data for NGC 4631 reveals strong, highly 
ordered, poloidal magnetic fields that extend several kiloparsecs above the galactic plane (Hummel et al 1991; GoUa 
& Hummel 1994; Beck 2000). These properties are by no means ubiquitous. In NGC 4565, the magnetic fields are 
strongly confined to the disk plane (Sukumar & Allen 1991). NGC 891 appears to be an intermediate case. Fields 
extend above the galactic plane, but they are not as highly ordered as those in NGC 4631 and the toroidal components 
are comparable to the poloidal ones (Hummel, Beck, & Dahlem 1991; Sukumar & Allen 1991). 

Recall that in SO configurations, the toroidal field dominates (the poloidal field vanishes in the symmetry plane). 
Thus the presence of a vertically averaged poloidal field that is strong relative to the azimuthal field may be indicative 
of an AO-type configuration. Sokoloff & Shukurov (1990) suggested that AO-type fields are generated by dynamos 
that operate in the gaseous halos of spiral galaxies. In disklike systems, there is a high gradient energy cost associated 
with AO-type field configurations. In spheroidal systems, this cost is reduced and therefore AO modes can grow and 
even dominate. Sokoloff & Shukurov (1990) proposed that two essentially independent dynamos operate in spiral 
galaxies, one in the disk and the other in the halo. Their order of magnitude estimates demonstrated that the model 
can explain the observations of NGC 4631 and NGC 891. Brandenburg et al. (1992) studied disk-halo dynamos in 
detail using numerical simulations. Their results confirmed, at a qualitative level, that a mean-field dynamo can 
operate in a turbulent, differentially rotating, spherical halo. One interesting property of the halo fields is that they 
are oscillatory (not unlike the solar cycle) with a period of oscillations comparable to the age of the galaxy. The ratio 
of odd to even parity fields therefore varies with time and is also sensitive to initial conditions. However, they found 
that it is generally difficult to construct models in which poloidal fields dominate to the extent seen in NGC 4631. 
Brandenburg et al. (1993) considered an alternative model in which outflows from the galactic disk due to galactic 
fountains, chimneys, and the Parker instability carry magnetic field from the disk to the halo. These fields are 
generally poloidal and might explain the magnetic structures observed in NGC 4631. They argued that by contrast. 
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the halo dynamo model works well for the galaxies NGC 891 and NGC 4565. 

As discussed above, spiral galaxies are found with axisymmetric, bisymmetric, and mixed symmetry field config- 
urations (Beck et al. 1996). In axisymmetric models, the m = mode is always preferred which suggests that some 
form of symmetry breaking is required in order to explain the observations. 

Even in axisymmetric models, growing mode solutions with m > are still possible. The generation of axisymmetric 
and bisymmetric fields was studied by Baryshnikova, Shukurov, Ruzmaikin, & Sokoloff (1987) and Krasheninnikova, 
Ruzmaikin, Sokoloff, & Shukurov (1989). They solved the three-dimensional dynamo equations for a thin disk by 
constructed a perturbation series where the small parameter is essentially the scale height of the disk divided by its 
scale radius. The growth rates for the two modes depend on the rotation curve one assumes for the disk and this 
result may explain the variations in field configurations observed in real galaxies. Results assuming a rotation curve 
appropriate to the galaxy M51 show that the characteristic growth timescale for the m — 1 bisymmetric mode can 
be as short as 2 x lO^yrs. This timescale is still a factor of four smaller than the growth rate for the axisymmetric 
mode. In addition, the field appears to be confined to a relatively narrow range in radius. 

Mestel and Subramanian (1991) considered a modified dynamo with a non-axisymmetric a-effect. Specifically, they 
imposed a (/)-dependence for the a-effect in the form of a uniformly rotating spiral structure. The motivation is the 
conjecture that the a-effect is enhanced along spiral arms. In the spiral density wave model a shock develops in 
the interstellar gas along the spiral arms. The jump in vorticity across the shock can lead to an enhancement of 
the a-effect which is constant in time but azimuth-dependent. As expected, an a-effect with a ^-dependence of this 
type naturally excites bisymmetric fields. The analysis of Mestel & Subramanian (1991) made use of the thin disk 
approximation and found rapidly growing bisymmetric magnetic fields that can extend over a sizable range in radius. 

H. Variations on the Dynamo Theme 

The dynamo hypothesis for the origin of galactic magnetic fields faces serious challenges beyond those raised above. 
In particular, observations have pushed back the epoch at which galactic-scale fields are known to exist thus shortening 
the time dynamo processes have to amplify small-amplitude seed fields. These difficulties have led to a number of 
variants of the standard aa;-dynamo. Here we discuss just a few. 

I. Parker Instability 

Parker (1992) proposed an alternative model that may alleviate the problems mentioned above while retaining the 
attractive features of the aw-dynamo. In his model, the combination of a hydromagnetic instability and reconnection 
of oppositely-oriented magnetic field lines (neutral point reconnection) provide the diffusion and a-effect necessary for 
the dynamo. The gas-magnetic field-cosmic ray system in galaxies is unstable (Parker 1979 and references therein). 
Consider a magnetic field that is initially purely toroidal. If any vertical waviness develops in the field, gas will slide 
downward toward the galactic plane making the outward bulges more buoyant. Relativistic cosmic rays produced in 
OB associations and supernovae can also force loops of magnetic field into the galactic halo. The result is a pattern 
of close-packed magnetic lobes perpendicular to the disk plane (Figure 11(a)). Magnetic reconnection can free these 
loops from the original toroidal field (Figure 11(b)). The tension along magnetic field lines is reduced thereby enabling 
differential rotation, the Coriolis effect, and further instabilities to act upon both the poloidal loops and toroidal field. 

In classical mean-field dynamo theory, turbulence is presumed to arise from nonmagnetic phenomena. Since backre- 
action of the field on the turbulence is assumed to be unimportant (at least until equipartition strengths are reached) 
properties of the turbulence, as expressed through the a and (3 tensors, are independent of the magnetic field. A 
dynamo based on the Parker instability represents a departure from this scenario in that the magnetic field (specif- 
ically, buoyancy of magnetic flux tubes) drives the a-effect. But while the underlying physics of Parker's model is 
quite different from that of the classical dynamo, the effective equations are essentially the same. The generation of 
poloidal field from toroidal field as well as the elimination of 'unwanted' flux can be described by a and /?-like terms 
in an equation for the large-scale held. Moreover the timescale for these processes is similar to what is obtained in 
the standard scenario. 

Moss, Shukurov & Sokoloff (1999) considered a mean-fleld disk dynamo in which the a-effect is driven by the Parker 
instability. In particular, they incorporated two modifications into the standard dynamo equation (Eq. (^7|)). First, 
a is taken to be an increasing function of |B|/_Boq where B^q is the equipartition field strength. Thus, the a-cffect 
becomes stronger as the field reaches equipartition, i.e., as the buoyancy of the flux tubes increases. Second, the 
rotational velocity is supplemented by a term corresponding to motions away from the disk, again due to buoyancy, 
i.e., V = rLo{r)cf) H- Vbz Vb also depends on |B|/i3oq. Solutions obtained from the full three-dimensional dynamo 
equations were found to resemble those obtained in conventional models in which the a-effect comes from cyclonic 
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turbulence. Some important differences did emerge. In particular, it was found that with a buoyancy- driven a-effect, 
the azimuthal field was no longer confined to a thin region near the equatorial plane but rather extended into the 
halo region. 

A variant of the Parker's model was considered in a series of papers by Hanasz & Lesch (1993, 1997, 1998). In 
particular, they studied the dynamics of a individual flux tubes under the influence of gravity, pressure from the 
ambient medium, and buoyancy due to both the magnetic field and cosmic rays pointing out that cosmic ray pressure 
inside a flux tube will make it buoyant even if the magnetic field is weak. In fact, with a weak field, magnetic tension, 
which tends to limit the Parker instability, will be unimportant. The hypothesis is then that star formation in an 
early phase of a galaxy leads to an excess of cosmic ray pressure over magnetic pressure. It is the enhanced buoyancy 
of the flux tubes due to the cosmic rays that drives and Parker-type fast dynamo, thereby converting cosmic ray 
energy to magnetic field energy. 



2. Magnetorotational Instability 

Magnetic fields may play a more direct role in the onset of turbulence. Balbus & Hawley (1991) have shown that a 
powerful, local, shear instability occurs in weakly magnetized differentially rotating disks. Almost any small seed field, 
in combination with an angular velocity profile that decreases with radius, will lead to dynamical instability. While 
the instability has been known for some time (Chandrasekhar 1960) it was Balbus & Hawley (1991) who recognized 
its importance for astrophysical disks (see Balbus & Hawley 1998 for a review). 

To understand the origin of the instability, consider an axisymmetric rotating disk with angular velocity lo — lij{R) 
and magnetic field parallel to the spin axis. Focus on two fluid elements along the same field line and consider what 
happens when they are displaced, one in the direction of rotation and the other in the opposite direction. Magnetic 
tension causes the leading element to slow and move inward and the trailing element to speed up and move outward. 
If there is differential rotation with du / dR < the fluid element that flows inward will enter a more rapidly rotating 
region of the disk and its angular velocity increases stretching the field line is stretched even further. This positive 
feedback loop is the essence of the instability (Balbus & Hawley 1991). 

The instability can act as one component of a magnetohydrodynamic dynamo (e.g., Brandenburg et al. 1995; 
Hawley, Gammie, & Balbus 1996). The flows due to the instability regenerate magnetic fields which in turn reinforce 
the turbulence. A surprising feature of this system is that the magnetic energy is somewhat higher than the kinetic 
energy of the gas, i.e., super-equipartition field strengths are achieved. 

It is not clear whether the Balbus-Hawley instability is relevant to galactic magnetic fields. The underlying cause 
of the instability comes from the tendency of a weak magnetic field to enforce corotation. Since the inner regions 
of a spiral galaxy are typically in approximate solid body rotation a necessary ingredient for the instability is all 
but absent. However, the Balbus-Hawley instability may be important in the creation of the magnetic fields that 
seed the dynamo. It is possible, for example, that seed fields are generated in active galactic nuclei (AGN) and then 
dispersed through the ISM by jets with the Balbus-Hawley instability playing a central role in the generation of these 
AGN-fields. 



3. Cross-Helicity Dynamo 

The cross-helicity dynamo is a variant of the standard mean-field dynamo which allows for more rapid amplification 
of a seed field at early times. The effect was first proposed in a general context by Yoshizawa (1990) and later applied 
to the problem of galactic magnetic fields by Yokoi (1996) and Brandenburg & Urpin (1998). The model assumes 
that a cross-correlation exists between the fluctuating components of the velocity and magnetic field, i.e., (v • b) ^ 0. 
The mean-field MHD equation is then supplemented by a term that is proportional to the product of (v • b) and the 
large-scale vorticity field, ^ = V x V. Since this term is not proportional to the large-scale magnetic field, it leads to 
linear rather than exponential growth of B. Under certain assumptions, the cross-helicity dynamo can dominate at 
early times speeding up the amplification process (Brandenburg and Urpin 1998). 

It is interesting to compare the cross-helicity effect with the standard a mechanism. In the cross-helicity dynamo, 
the mean vorticity distorts V through the inertial term in the Euler equation (Eq. (p^)) 



_ = _(vV)V + ... , (80) 
and distorts b through the stretching term in the MHD equation (Eq. (||)) 
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By combining Eqs. (pQ) and (p^) we find 



^ = (b-V)V+... . (81) 



dt 



(v X b) = b X (v V) V + V X (b • V) V - 



(82) 



If the fluctuation fields are isotropic, (vibj) ~ 
force of tlie form Sen = XC+ = ■ ■ ■ wliere A ; 



i(5y (v • b). Tlie cross- Irelicity effect tlrerefore leads to an electromotive 
E |r(v • b) and r is the characteristic correlation or turnover time of the 
turbulence. When this term is included, the mean-field dynamo equation takes the form 



— = V X (V X B) + V X (aB - /3V X B + ^ch) ■ (83) 
at ' ' 

Thus, there is a source term S = V x £ch in the induction equation which is independent of the large-scale magnetic 
field. The effective current associated with this source is parallel to the mean vorticity which, in a disk galaxy, is 
parallel to the spin axis. Only toroidal fields are generated. By contrast, the a-effect generates a current which is 
parallel to the magnetic field and hence poloidal fields are generated from toroidal ones. 

Brandenburg & Urpin (1998) considered the cross-helicity dynamo in the context of galactic magnetic fields. In 
order to approximate the importance of the effect, they analyzed different data sets from MHD turbulence simulations 
and concluded that the relative cross-correlation parameter, e = (v • b)/ ((v)(b)) is in the range 3 x 10~^ — 3 x 10~^. 
Using values characteristic of spiral galaxies S can be parametrized as follows: 




The cross-helicity dynamo dominates over the standard aw-dynamo for values of the field strength B < ST ^ where, 
as before, F is the growth rate of the fastest growing dynamo mode. For example, if F ~ 0.5 Gyr^^, the transition 
from cross-helicity to aw-dynamo occurs at a field strength i? ~ 5 x 10~® — 5 x 10~^° (depending on the value used 
for e) and at a time t ~ 1 Gyr). The cross-helicity effect is potentially important in cases where the initial seed field is 
extremely weak. In the standard dynamo model, a seed field of strength 10~^* G is amplified to a strength of 10~* G 
in ~ 4.5 Gyr. In a cross- helicity/aw-hybrid dynamo, the same field strength can be reached in ~ 1 Gyr (Brandenburg 
& Urpin 1998). 



I. Dynamos in Irregular and Elliptical Galaxies and Galaxy Clusters 

The existence of microgauss fields in elliptical galaxies and galaxy clusters presents distinct challenges to models for 
the origin of cosmic magnetic fields. Rotation plays a central role in the mean-field dynamo models devised to explain 
magnetic fields in disk galaxies, stars, and planets. Rotation provides a resevoir of energy for field amplification 
through both a and uj effects. In particular, the a-effect requires net helicity, or equivalently, mirror symmetry 
violation in the turbulence. In rapidly rotating systems mirror symmetry violation occurs because of the Coriolis 
effect. This effect is weak, if not absent, in slow rotating systems such as ellipticals and clusters. 

A turbulent dynamo may well operate in ellipticals and clusters, though without rapid rotation the coherence length 
of the fields is limited by the characteristic scale of the turbulence. In this section, we discuss the necessity of dynamo 
action in ellipticals and clusters, sources of energy for the generation of turbulence, and the interplay between cooling 
flows and magnetic flelds. 
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1. Elliptical Galaxies 

Moss & Shukurov (1996) identified two potential sources for seed fields in elliptical galaxies: stellar magnetic fields 
ejected into the ISM by supernavae and stellar winds and magnetic remnants that arise if ellipticals formed from 
mergers of spiral galaxies. The former can lead to fields of about 10~^^ G (see below) while the latter can lead 
to fields or about 10~^ G. In either case, Moss & Shukurov (1996) concluded that the observed microgauss fields in 
ellipticals require further amplification. Moreover, simple stretching of fields lines by plasma motions without dynamo 
action is probably unrealistic since Ohmic dissipation will lead to loss of field energy. They therefore concluded that 
magnetic fields arc amplified by dynamo processes. However, the conditions in elliptical galaxies do not appear to 
support a mean-field dynamo for the reasons discussed above. (See, however, Lesch & Bender (1990) where it is 
argued that a mean-field dynamo can operate in elliptical galaxies.) Instead, Moss & Shukurov (1996) proposed that 
fields in ellipticals arc amplified by a so-called fluctuation dynamo (Kazantsev, Ruzmaikin, & Sokoloff 1985) in which 
turbulent motions lead to rapid growth of the rms magnetic field. 

Moss & Shukurov (1996) identified on two potential sources of turbulence in ellipticals, ejecta from supernovae and 
stellar winds and the stirring of the interstellar gas by the random motions of stars. The former gives rise to acoustic 
turbulence — essentially the random superposition of sound waves — with a characteristic scale ~ 1 kpc. The latter 
can drive vortical turbulence with a much smaller characteristic scale (~ 3pc). Moss & Shukurov (1996) suggested 
that amplification of magnetic fields in ellipticals takes place in two stages. First, vortical turbulence, acting as a 
dynamo, amplifies a seed field to cquipartition strengths on 2 — 3pc scales within a very short (< 10'^ yr) timescale. 
Acoustic turbulence amplifies the field and also leads to an increase in the coherence length. The model predicts 
that the fields in the inner regions of ellipticals will have strengths that range from 14 /xG near the center to 1 /xG at 
r = 10 kpc and have a characteristic scale of 100 — 200 pc. 

The effects of cooling flows as well as a time-dependent supernovae rate were considered by Mathews & Brighenti 
(1997). Moreover, they found that compression of the magnetic field in cooling flows could lead to an additional 
ampliflcation by a factor of 10^. Moreover, the exact value of the field in ellipticals is sensitive to interstellar turbulence 
in the distant past. 

2. Clusters 

The origin of magnetic fields in clusters is, at present, uncertain. Cluster galaxies are an obvious source of magnetic 
fields. Material originally associated with individual galaxies is spread throughout the intracluster medium by tidal 
stripping and galactic outflows. Magnetic fields, tied to the fluid, are likewise dispersed throughout the cluster, 
albeit with a smaller field strength due to the usual V^^'^ dilution factor. The observation that fields in clusters are 
comparable in strength to those in galaxies implies that field amplification is occurring in the intracluster medium. 

Early attempts to understand the amplification of magnetic fields in clusters focused on dynamo action in galactic 
wakes (Jaffe 1980; Roland 1981; Ruzmaikin, Sokoloff, & Shukurov 1989). As galaxies move through the intracluster 
medium they generate a turbulent wake which can presumably support a dynamo. However, more detailed calculations 
have shown that it is very difficult to produce fields above ^ 10~^ G by this mechanism (Goldman & Rephaeli (1991) 
and De Young (1992)). 

A more promising scenario relies on mergers to drive cluster-scale dynamos (Tribble 1993). Clusters are relatively 

young systems and most have undergone at least one merger event during the last Hubble time. The amount of 
energy released during a major merger is comparable to the gravitational and thermal energies of the system and 
larger than the magnetic energy by a factor ^ 10^ — 10'^. Therefore only a small fraction of the available energy is 
required to explain cluster fields. Mergers lead naturally to turbulence and shocks in the intraclustc^r medium. The 
former is necessary for a dynamo while the latter can accelerate particles which can then produce the synchrotron 
emission. The lifetime of magnetic fields in the intracluster medium can be relatively long (Soker & Sarazin 1990; 
Tribble 1993). This timescale is set by the rate of magnetic reconnection, which in turn is limited by the Alfven speed. 
For a B ~ 1 /zG field and electron number density rie ~ 2 x 10~^cm~^ one finds an Alfven speed va — 30 kms~^ so 
that the reconnection time for 10 kpc magnetic structures is > 3 x 10^ yr. (Note, however that the timescale for the 
field to decay is shorter if the field is concentrated in thin flux tubes (Subramanian 1998; Ruzmaikin et al. 1989).) By 
contrast, the lifetime of synchrotron emitting electrons is short. Relativistic electrons lose energy by inverse Compton 
scattering off CMB photons. The lifetime at the present epoch is ^ 10^ yr and decreases rapidly with redshift (Tribble 
1993). Thus, radio halos fade quickly while cluster magnetic fields survive for periods comparable to the Hubble time. 
This simple conclusion may explain the rarity of radio halos and the ubiquity of cluster magnetic flelds (Tribble 1993). 

As in elliptical galaxies, cooling flows can have a significant effect on the magnetic fields in clusters. Magnetic 
fields in cooling fiows are amplified by radial infall and shear (Soker & Sarazin 1990) and can be expected to reach 
equipartition with the gas. A direct observable effect of these fields is likely to be very strong Faraday rotation that 
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increases rapidly toward the cluster center. 
V. SEED FIELDS 

Though magnetic fields could have been a feature of the initial conditions of the Universe a more appealing hypoth- 
esis is that they are created by physical processes operating after the Big Bang. These first fields can be extremely 
small since subsequent dynamo action can amplify them by many orders of magnitude. However, even small fields 
require an explanation (Zcl'dovich & Novikov 1983; Rces 1987; Kronberg 1994; Beck et al. 1996; Kulsrud 1999). 

The list of proposals for the origin of seed fields is now long and diverse. However no single compelling model has 
emerged and thus the following questions remain: 

• When did magnetic fields first appear? Were they present during big bang nucleosynthesis? recombination? 
galaxy formation? 

• What was the spectrum (strength vs. coherence length) of the first magnetic fields? Was the galactic dynamo 
seeded by subgalactic, galactic, or supergalactic-scale fields? In the language of the old debate on cosmological 
structure formation, are galactic magnetic fields a top-down or bottom-up phenomena? 

• Is there a connection between the creation of the first fields and the formation of large-scale structure? 

• Evidently, a dynamo is necessary to maintian galactic magnetic fields. Is dynamo action also necessary for 
amplification of an initially small seed field? 

Scenarios for seed fields fall into two broad categories, those that rely on ordinary astrophysical processes and 
those that rely on new and exotic physics. By and large, astrophysical mechanisms exploit the difference in mobility 
between electrons and ions. This difi'erence can lead to electric currents and hence magnetic fields. Mechanisms of 
this type can operate during galaxy formation or alternatively they can operate in other systems such as stars and 
active galactic nuclei. In the latter case, the question is to explain how magnetic fields in those systems make their 
way into the interstellar or protogalactic medium. 

Exotic processes in the very early Universe can also create magnetic fields. There is strong circumstantial evidence 
that the Universe has undergone a series of phase transitions since the Big Bang including an episode of inflation and 
the electroweak and quark- hadron transitions. Magnetic fields of interesting strength can arise during these events 
though their connection to galactic and extragalactic magnetic fields remains unclear. 

A. Minimum Seed Field for the Galactic Dynamo 

What is the minimum seed field required to explain the observed galactic and extragalactic magnetic fields? Un- 
certainties in our understanding of the dynamo, galaxy formation, and cosmological parameters make this quantity 
difficult to pin down. In general, the dynamo amplification factor A can be written 

A=^ = e^^*/-**) (86) 

where F is the growth rate for the dominant mode of the dynamo, Bi is the field strength at the time U when the 
dynamo begins to operate, and i/ is the time when the fields reach the observed value Bf. 

Irrespective of dynamo action, a seed field that is created prior to galaxy formation is amplified as the protogalactic 
gas collapses to form a disk. This point was stressed by Lesch & Chiba (1995) who estimated the amount of pre- 
dynamo amplification that occurs during galaxy formation. Their analysis was based on a simple and well-worn model 
for the formation of a disk galaxy proposed by Fall & Efstathiou (1980) (see, also White & Rees 1978). According 
to Fall & Efstathiou (1980), the first stage of galaxy formation is the development of an extended virialized halo of 
gas and dark matter where the gas and dark matter are assumed to have the same specific angular momentum. Gas 
elements cool and lose energy but, by assumption, conserve angular momentum. The end result is a disk which, for 
simple and reasonable choices of initial conditions (e.g., a gas-dark matter halo in solid-body rotation), resembles the 
disks of present-day galaxies. 

In the hierarchical clustering scenario, halos form around peaks in the primordial fluctuation distribution. Small- 
scale objects collapse first and coalesce to form systems of increasing size. In general, the inner regions of halos 
form first. The spherical infall model (Gunn & Gott 1972) provides a useful, albeit highly idealized, picture of halo 
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formation. In this model, the mass distribution around a given peak is assumed to be spherically symmetric with a 
density profile that decreases monotonically with distance from the peak center. The evolution of a protohalo can 
therefore be described in terms of the dynamics of spherical shells. Initially, the shells expand with the Hubble flow, 
though eventually they reach a maximum radius known as the turnaround radius Rta, break away from the general 
expansion and collapse. After a shell collapses, it virializes with the other shells (though this process is not particularly 
well understood). The model captures certain features of the hierarchical clustering scenario (e.g., inside-out halo 
formation) while ignoring the complexities of true hierarchical structure formation. A number of results obtained 
from this model (e.g., characteristic density, halo formation time) are in good agreement with those obtained from 
N-body simulations of hierarchical structure formations. Simple arguments suggest that the characteristic radius of a 
virialized shell is roughly one half of the turnaround radius, i?vir — -Rta/2- For a typical spiral galaxy, i?vir — 150 kpc. 
A seed magnetic field created prior to halo formation is therefore amplified by a factor (i?ta/^vir)^ — 4. 

A crucial assumption for what follows is that the gas-dark matter halo acquire angular momentum as it forms. In 
a hierarchical clustering scenario, protogalaxies gain angular momentum through tidal interactions with neighboring 
protogalaxies (Hoyle 1949; Peebles 1969; White 1984). (Tidal torque theory, of course, requires a departure from 
spherical symmetry since tidal fields do not couple to spherical shells.) It is common practice to characterize the 
angular momentmn of a system in terms of the dimcnsionlcss spin parameter A = J\E\^/^ /GM^^^^ where J is the total 
angular momentum of the system and E is its total energy. Roughly speaking, A is the square root of the ratio of 
the centripetal acceleration to the gravitational acceleration. Theoretical analysis and numerical simulations suggest 
that, for a typical halo, A ~ 0.07 (Peebles 1969; Barnes & Efstathiou 1987). 

As described above the second stage in the formation of a disk galaxy, the gas dissipates energy and sinks to the; 
center of the halo's gravitational potential well where it forms a rotationally supported disk (i.e., a system in which 
A ~ 1). If the protogalaxy were composed entirely of gas, the radius of the disk would be i?disk — A^i?vir ^ 1 kpc. 
Not only is this value too small (the typical scale radius of a spiral galaxy is ^ 10 kpc) but the time required 
for such a the disk to form is greater than the age of the Universe (Peebles 1993). In a CDM model, 90% of a 
protogalaxy is in the form of collisionless dark matter which remains in an extended halo. In this case, the disk radius 
is i?disk — Ai?vir — 10 kpc. 

Magnetic fields, frozen into the gas, are amplified by a factor A~^. Additional amplification occurs due to the 
collapse along the spin axis as the gas forms a thin disk. The net pre-dynamo amplification factor is therefore 



^' ~ 8 X 10^ ^J^taARvir^ ^£i^iy ^ Rdisk/ Rh 



(87) 



where Rh is the disk scale height (Lesch & Chiba 1995). 

The quantity T is model-dependent and not very well known. In Section IV. E, we found 



where c = 0.6 — 0.8 depending on the functional form of a{z). Dimensional arguments suggest that = Z/^w/Zi where 
L is the size of the largest eddies and h is the scale height of the disk (Ruzmaikin, Shukurov, and Sokoloff 1988a). 
With h ~ 500 pc, w ~ 29kms~ikpc, and L ~ 100 pc, we find ao - 6 x lO-^cms'^ and T < 0.7- 0.9Gyr"^ Ferriere 
(1992, 1993) derives the somewhat smaller value a ~ 2 x lO^cms"^ for the supernovae-driven model of turbulence 
and detailed models of a galactic dynamo based on these results give F ~ 0.45Gyr~^ (Ferriere and Schmidt 2000). 
However, as discussed above, their calculation may underestimate the strength of the a-effect from supernovae. 

Observations of CMB angular anisotropy, high-redshift supernovae, and large-scale structure indicate that the 
Universe is spatially fiat with flm — 0.15 — 0.40 and J^a — 1 — flm- This conclusion implies an older Universe as 
compared with one in which f2„ = 1. The lower bound on the required seed field is therefore relaxed since there is 
more time for the dynamo to operate. In a spatially flat, Robertson- Walker Universe, 
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where, as before, a = (1 + z)~^ is the scale factor at time t with a(io) = 1 and A = {Q\/il.m)^^^ (1 + z)~'^ . An upper 
bound on Zi is given by the redshift at which a protogalaxy separates from the Hubble flow, collapses, and virializes. 
Simple estimates based on the spherical infall model suggest a value Zi ~ 50. However, disks are almost certainly 
assembled at a much later epoch. Observerations such as the Hubble Space Telescope Deep Field Survey, for example, 
indicate that disk galaxies at z ~ 3 are highly irregular and still in the process of being formed (see, for example, 
Lowenthal et al. 1997) 

A plot of tf — ti vs. for Hq = 70kms~^Mpc, tf = t^, and various choices of Zi is given in Figure 12(a). The 
amplification factor A, assuming a growth rate F = 2Gyr~^, is shown on the right-hand vertical axis. Note that 
log A scales linearly with T/Hq. Evidently, A varies by several orders of magnitude depending on Zi and ^Im- For 
Zi = 10 and flm = 0.2, the maximum amplification factor is ^ ~ 10^^ and therefore a present-day microgauss requires 
a seed field with strength Bi ~ 10^^° G. (We have not included the amplification A' which can occur during galaxy 
formation.) 

Davis, Lilley, & Tornquist (1999) argued that the lower bound on the strength of the seed field is 10^'^*' G or 
less depending on the cosmological model. Their estimate assumes a lower value for the Hubble constant {Hq = 
50 — 65kms~^Mpc), values for the dynamo growth rate as high F — 3.3Gyr~^, and an earlier choice for ti. The 
latter value is taken to be the time when a region that corresponds in size to the largest turbulent eddy 100 pc) 
breaks away from the Hubble flow and collapses. The implicit assumption is that dynamo processes begin to operate 
on subgalactic scales well before the disk forms. While this assumption is reasonable, it requires further study. In 
particular, the estimates for F, derived for a disk dynamo, do not necessarily apply to subgalactic objects in an 
evolving hierarchy of structure. 

Observations of microgauss fields in galaxies at moderate redshifts tighten severely the lower bound on seed fields. 
The results of Kronberg, Perry, and Zukowski (1992) imply Zf > 0.4 while those of Wolfe, Lanzetta, and Oren (1992) 
suggest a more stringent albeit controversial choice Zf > 2. zj >2 is also implied by the observations of radio galaxies 
by Athreya et al. (1998). Figures 12(b) and 12(c) show tf — ti for Zf = 0.4 and 2 respectively. With Zf = 0.4, Zi = 10, 
and = 0.2, the limit on Bi is tightened to 10~^^G. More severe is the case Zf = 2 where the time available for 
the dynamo to operate is shortened to a few billion years or less and a seed field with Bi ~ 10^ ^'^ G is required. 



B. Astrophysical Mechanisms 

The difference in mobility between electrons and ions in an ionized plasma leads to charge separation effects and a 
breakdown of the MHD approximation. Gonsider a multicomponent fluid composed of electrons, protons, and photons. 
(The discussion can be extended easily to include heavier ions and neutral atoms). The momentum equations for the 
electrons and protons can be written: 



_ E -f ^^-^ + (90) 



dt Pe Top V C I TO 



dVp _ Vpp e ( Yp X B\ Kep 



dt 



Pp 



TO„ 



E + _ V(/. - =^ (91) 



where pi, rrii, pi = mim, Ui, and are the partial pressure, mass, mass density, number density, and velocity of 
the electrons {i — e) and protons {i ~ p) (e.g., Spitzer 1962; Sturrock 1994). The rate of momentum transfer from 
protons to electrons is Kgp = TOe (Vp — Vg) /r where r = TOeCr/nge^ is the characteristic timescale of electron-proton 
collisions. The current density is J = e (upVp — rigVe) and therefore, if approximate local charge neutrality holds 
(i.e., He ^ rip), Kep ~ eJ/cr. 

Also important for some scenarios is the coupling of electrons and photons due to Thomson scattering. In the strong 
coupling limit, photons behave like an ideal fluid. Thomson scattering leads to an additional momentum transfer term 
in Eq. ( |90| ) of the form K.^-/ = '^caxp-y (V^ — Ve) /3 where ctt = fTJ" (e^/TOgC^)^ = 6.65 x 10~^^ cm^ is the Thomson 
cross section and p-yC^, p-y, and V-y are respectively the energy density, pressure, and (fluid) velocity of the photon 
fluid. The momentum equation for the photon fluid is given by 
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(see, for example, Peebles 1980). 
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«^K. + :^B ^— — a^V^B (93) 



1. Seed Fields from Radiation-Era Vorticity 

An early attempt to explain the origin of seed fields is due to Harrison (1970, 1973) who considered the evolution 
of a rotating protogalaxy prior to decoupling. During this epoch, electrons and photons are tightly coupled so that, 
for most purposes, they can be treated as a single fluid. The coupling between electrons and protons is somewhat 
weaker allowing currents and hence magnetic fields to develop. 

Harrison considered a homogeneous, spherical protogalaxy in solid body rotation with radius R. As the protogalaxy 
expands, PpR^ and p-yR'^ remain constant. In the limit that the electron-proton coupling is ignored the angular 
momenta of the electron-photon and proton fiuids are separately conserved, i.e., ppUj^R^ and p^us'^R^ are constant. 
This implies that ujp oc R~^ and uj^ oc R^^. The ions "spin down" faster than the electrons and photons and a current 
J ~ euj^R/nfipC and hence a magnetic field B ~ mpCUJ^/e ~ 10~^ G (w-y/s"^) develop. 

A more formal derivation of this result, including the electron-proton coupling, follows from Eqs. (po|)-(|9^). The 
velocity field for Harrison's spherical protogalaxy can be decomposed into a homogeneous expansion and a circulation 
flow, i.e., Vp = (a/a) r -I- Up where (^^ = V x Vp = V x Up is the vorticity of the proton fluid. Combining the curl of 
Eq. (|9^) with Maxwell's equations gives 

d 

dt \ ttihc J J ATTarriH 

where, because of the assumed homogeneity, the pressure gradient term in Eq. (|9^) is dropped. In the highly conducting 
protogalactic medium, the diffusion term is negligible thus yielding the conservation law 

C ^t) + ^Bit)=(^)\p{U) (94) 
niHC \ a[t) J 

where ti is some initial time characterized by zero magnetic field. 

The dynamics of the electrons is determined primarily by the photons. The inertial and gravitational terms in 
Eq.(pO[) are therefore neglected leaving what is essentially a constraint equation for the electron velocity field: 

E+Xl^ = £^ + ^^. (95) 
c a e 

The curl of this equation, together with Maxwell's equations, gives 

1^ (a^B) = X Ke^ (96) 

a'^ dt ^ e 

while the curl of Eq. (63) gives 



a dt ^"^^^ 4p^7 

Eqs. (|9^), (|96|), and ( |97| ) can be combined to eliminate K.^^ and Cp(^i) in favor of B 
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In deriving this equation we have used the fact that p^ cx a^^ and pp oc a~^. In addition, to leading order, we set 
= Ve = Vp and C-y = Cpj i-S-, differences between the electron, proton, and photon velocities enter the calculation 
only through the collision terms. The desired result, B ~ nipcQ/ e, follows immediately if we take a{t) ^ a{ti). By way 
of example, we note that the solar-neighborhood value C, ~ 29kms~^ kpc~^ yields a seed magnetic field B ~ 10^^^ G. 

Harrison's scenario for the origin of magnetic fields has a number of attractive features which are echoed in subse- 
quent proposals. First, seed fields are generated as part of the galaxy formation process and therefore naturally have 
a coherence length comparable to the size of the galaxy. Second, the seed fields have a strength that is set by the 
vorticity of the protogalaxy which can be related to the present-day vorticity in galactic disks, an observable quantity. 
These fields, though weak, are sufficient to explain present-day galactic fields, provided an efficient dynamo develops 
to amplify them. The most severe criticism of the model is that, prior to structure formation, vorticity decays rapidly 
due to the expansion of the Universe (Rees 1987). The implication is that vorticity in disk galaxies in not primordial 
but rather generated during structure formation. Harrison's scenario does suggest that there may be a connection 
between the origin of galactic vorticity and the origin of seed fields. 



45 



2. Biermann Battery Effect 

In the hierarchical clustering scenario, protogalaxies acquire rotational angular momentum from the tidal torques 
produced by their neighbors (Hoyle 1949; Peebles 1969; White 1984). However, gravitational forces alone cannot 
generate vorticity and therefore its existence in galactic disks must be due to gasdynamical processes such as those 
that occur in oblique shocks. In an ionized plasma, these same processes produce magnetic fields (Pudritz & Silk 
1989; Kulsrud, Cen, Ostriker, & Ryu 1997; Davies & Widrow 2000). 

We consider vorticity generation first because it is conceptually simpler. In the absence of electromagnetic effects 
and viscosity, the evolution of a coUisional fluid is given by Eq. (jlj) with J = and j] ~ 0. Taking the curl yields the 
following equation for the vorticity: 



|-VX(VXC, .Ii!2Ll£. (99) 

The source term on the right-hand side comes from gasdynamical effects, namely pressure and density gradients that 
are not coUinear. (The gravitational force term does not appear since the curl of VV^ is identically zero.) Vorticity 
generation is illustrated in Figure 13 where we consider an ideal single component fluid consisting of a particles of 
mass m. The gravitational force mg is only partially compensated by the pressure gradient force — Vp/n. The 
acceleration rate for high density regions is greater than the rate from low density regions and therefore the velocity 
field downstream of the pressure gradient has a nonzero vorticity and shear [dvy/dx ^ 0). 

In an ionized plasma, similar inertial effects lead to electric currents and magnetic fields. This mechanism, known 
as the Biermann battery effect (Biermann 1950; Roxburgh 1966) can be derived by combining the generalized Ohm's 
law with Maxwell's equations. The former, found by taking the difference of Eqs. (|9^) and ( pi|) can be written 
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where terms of order me/mp^ as well as terms quadratic in the velocities, have been neglected (Spitzer 1962; Sturrock 
1994). In addition, we have assumed that local charge neutrality (rie = Up) holds and that the electron and proton 
partial pressures are equal. The J x B term describes the backreaction of the magnetic field on the fluid. Since, 
in the present discussion, we are interested in the creation of (small) seed fields, this term can be ignored. In ideal 
MHD, the left-hand side, which described charge separation effects, is zero. In addition, the pressure gradient term 
is ignored. We therefore recover the simple form of Ohm's law: J = cr (E -I- (Vp x B) /c) where Vp is essentially the 
fluid velocity. In the present discussion, we focus on charge separation effects and assume that the conductivity is 
high and that B and Vp/c are small. The result is a form of Ohm's law which is relevant to the Biermann effect: 
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The curl of this equation. 
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together with Maxwell's equations, gives 



gg-Vx(VxB)^^ ^^-"7^- . (103) 
at e 

The Biermann effect can be illustrated by the following simple example. (For a similar pedagogical discussion 
of the Biermann effect in stars, see see Kemp (1982)). Consider, first, the plane-symmetric flow in Figure 14(a). 
Since the gravitational force on the protons is much greater than that on the electrons, electrostatic equilibrium 
{dVp/dt = dVe/dt) requires an electric field E ~ —^p/eue — mg'/2 where g' = g — dVp/dt is the effective 
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gravitational force in the rest frame of the fluid. An electrostatic force of this type exists in stars (Rosseland 1924) 
but is extremely small and inconsequential to understanding stellar structure. 

Electrostatic equilibrium is no longer possible when gradients in theremodynamic quantities such as the density and 
temperature are not parallel to the pressure gradient. This situation is illustrated in Figure 14(b) where the density 



increases in the +?/-dircction. The second term in Eq. (102) is a nonzero vector in the +z-direction and therefore at 



least one of the other terms in this equation is nonzero. In fact, since a nonzero V x E implies a time-dependent 
magnetic field, which in turn implies time-dependent currents, both terms are nonzero. However, the term on the 
left-hand side is negligible: 
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where L i s th e typical length scale for the system. We therefore find a time-dependent magnetic field, in agreement 



with Eq. (103). For the situation shown in Figure 14(b), the magnetic field is in the -t-z-direction. 

Approximate local charge neutrality implies that rig ~ rip = xp/rup where x is the ionization fraction. In addition, 
since the electrons and protons are expected to be in approximate thermal equilibrium, pe cii prie/ {ric + rip) — 
Pxl (1 + x) where p is the total gas pressure. The source term on the right-hand side of Eq. ( |103| ) can therefore be 
written 
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and is evidently proportional to the source term for the vorticity. 

It is interesting to note that the Biermann effect has been observed in lasar-generated plasmas (Stamper & Ripin 
1975; also see Loeb & Eliezer 1986). The schematic diagram of a typical experiment is shown in Figure 15. The 
plasma has strong temperature and density gradients that are nearly perpendicular leading to a source term for the 
magnetic field. These magnetic fields, typically megaguass in strength, were studied through Faraday rotation induced 
on a second laser beam that was shined through the plasma. 

Kulsrud et al. (1997) simulated the creation of magnetic fields via the Biermann battery effect. They used a 
cosmological hydrodynamic code that was designed to handle shocks as well as gravitational collapse (Ryu et al. 1993). 
Since the fields are extremely weak, their backreaction on the fluid could be ignored. The MHD equation, supplemented 



by the Biermann term (Eq. (103)) was solved with V obtained from the hydrodynamic code. Kulsrud et al. (1997) 
assumed parameters appropriate to the CDM model that was popular at the time of their work {h = 0.5, = 0.06 
and r^o = 1) and found that galactic-scale fields are created with strengths of order 10~^^ G. However, the spatial 
resolution of their simulations is ~ Mpc, comparable to, if not somewhat larger than, the scales of interest. 

Davies & Widrow (2000) took a complementary approach by focusing on the collapse of an isolated galaxy-sized 
object in an otherwise homogeneous and isotropic universe. Their investigation is in the spirit of semianalytic models 
of disk galaxy formation by Mestel (1963), Fall & Efstathiou (1980), and Dalcanton, Spergel, & Summers (1997) 
as well as numerical simulations by Katz & Gunn (1991). However, those models assume that angular momentum 
and vorticity are present ab initio and therefore they are unable to shed light on voriticity and hence magnetic field 
generation. In Davies & Widrow (2000) vorticity generation is followed explicitly. The system considered consists 
of coUisionless dark matter and coUisional gas with an initial perturbation whose density distribution is assumed to 
be nearly spherical, smooth (i.e., no small-scale perturbations) and monotonic, specifically, that of an axisymmetric, 
nonrotating prolate protogalaxy. As discussed above, each fluid element expands to a maximum or turnaround radius 
before falling in toward the center of the protogalaxy. During the early stages of collapse, an outward moving shock 
develops and as the infalling gas crosses the shock, it is heated rapidly and decelerated (Bertschinger 1985). In 
addition, the velocity changes direction at the shock. This is where both magnetic fields and vorticity are generated. 

Analytic calculations based on the thin-shock approximation (Landau & Lifshitz 1987) and numerical simulations 
using SPH follow explicitly the generation of vorticity at the shock that forms. In this simple model, the vorticity 
and hence magnetic field are in the azimuthal direction and are antisymmetric about the equatorial plane so that the 
total vorticity of the system is zero. This of course refiects that fact that angular momentum has not been included. 
Rotation about a short axis of the protogalaxy will "shear" the vorticity and magnetic field lines into a dipole-like 
configuration. 

The results of this analysis indicate that by a redshift z ~ 8, a lO^'^" G field, coherent on 10 — 20kpc scales, is 
generated. This field is amplified by a factor 10^ — 10'^ as the protogalaxy collapses (Lesch & Chiba 1995) leading to 
a 10~^^ G field in the fully formed disk galaxy. 
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In a hierarchical scenario, the Biermann effect also operates in subgalactic objects suggesting an alternate route to 
galactic-scale seed fields. Consider a region that contains a mass M ~ 10® Mq in baryons. This mass corresponds 
to the Jeans mass at decoupling and represents a lower bound on the mass of the first generation of gravitationally 
bound objects. Moreover, gas clouds of this size are thought to be the first sites of star-formation. The Biermann 
effect leads to seed fields of order IQ-^^ G by z ~ 40 (Pudritz & Silk 1989; Davies & Widrow 2000). The dynamical 
time for these objects is relatively short (< 10^ yr) and therefore amplification by a dynamo can be very fast. The 
field reach equipartion (microgauss) strength on these subgalactic scales. As discussed below, small-scales fields of 
this type can act as a seed for a galactic-scale dynamo. Moreover, the existence of magnetic fields in star-forming 
clouds may resolve the angular momentum problem for the first generation of stars (Pudritz & Silk 1987). 

Lazarian (1992) considered a variant of the Biermann battery in which electron diffusion plays the key role in estab- 
lishing electric currents and hence magnetic fields. This effect relies on the observation that the ISM is nonuniform, 
multiphase, and clumpy. Diffusion of electrons from high to low density regions results in an electric field whose curl 
is given by 



Lazarian (1992) estimated that magnetic fields as high as 3 x 10~^^ G can be generated on large scales in the ISM. 

A battery of this type may have operated prior to the epoch of galaxy formation (Subramanian, Narasimha, & Chitre 
1994). At decoupling [zd — 1100), protons and electrons combine to form neutral hydrogen. However, the absence 
of Lya absorption in quasar spectra (above and beyond the features in the Lya forest) severely limit the amount of 
smoothly distributed HI in the intergalactic medium. The implication is that hydrogen has been reionized almost 
completely at some epoch between decoupling and z ~ 5 (Gunn & Peterson 1965). Subramanian, Narasimha, & Chitre 
(1994) assumed that reionization is characterized by ionization fronts propagating through the intergalactic medium. 
Electric currents and magnetic fields are produced when these fronts encounter density inhomogeneities. Gnedin, 
Fcrrara, & Zwcibcl (2000) performed numerical simulations detailing this mechanism. These simulations followed the 
reionization of the Universe by stars in protogalaxies. Ionization fronts formed in high density protogalaxies propagate 
through the Universe. As they cross filamentary structures in the so-called cosmic web, they drive currents which in 
turn give rise to magnetic fields. The fields produced are highly ordered on a megaparsec scale and have a strength 
of 10^^^ G. Furthermore, the fields are stronger inside dense protogalaxies though a detailed analysis on these scales 
was limited by resolution of the simulations. 

3. Galactic Magnetic Fields from Stars 

The first generation of stars can be a source of galactic-scale seed fields. Even if a star is born without a magnetic 
field, a Biermann battery will generate a weak field which can be amplified rapidly by a stellar dynamo. If the star 
explodes or undergoes significant mass loss, magnetized material will be spread throughout the ISM. The following 
simple argument illustrates why this proposal is attractive (Syrovatskii 1970). Over the lifetime of the Galaxy, there 
have been ^ 3 x 10* supernova events, roughly one for each (10 pc)^ volume element. The magnetic field in the 
Crab nebula, taken as the prototypical supernova remnant, is 300 /iG over a region 1 pc in size. A galaxy filled with 
(ancient) Crab-like nebulae would therefore have an average field with strength ~ 3 jiG. Similarly, gas lost from stars 
via a wind will carry stellar magnetic fields into the ISM (Michel & Yahil 1973). 

The magnetic fields produced by stars will have a tangled component which is orders of magnitude larger than the 
coherent component. By contrast, the random and regular components of a typical spiral galaxy are nearly equal. 
Therefore, stellar magnetic fields, in and of themselves, cannot explain magnetic fields in disk galaxies. Nevertheless, 
the (rms) large-scale component due to an ensemble of small-scale stellar fields can act as a seed for a galactic dynamo. 

A naive model for the large-scale structure of stellar-produced magnetic fields divides the galactic disk into a large 
number of random cells of size i?ceii- The magnetic field is assumed to be uniform within each cell but uncorrelated 
from one cell to the next. The rms flux through a surface of scale i is F ~ BceiiiV^/^ ~ Bceii {L/Rce\\) where N is 
the number of cells that intersect the surface. The rms field on the scale L is therefore 
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where B^en is the average field in an individual cell. 

A more careful treatment takes into account the expected topology of the magnetic field. In particular, magnetic 
field lines within an individual cell should close. This point was stressed by Hogan (1983) in his discussion of the 
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magnetic fields produced in a cosniological phase transitions (see, also Ruzmaikin, SokolofFfe Shukurov 1988a, 1988b). 
Consider an idealized model in which the field within a given cell has a dipole structure. The orientation and strength 
of the dipoles are assumed to be uncorrelated from one cell to the next. Let us estimate the magnetic flux through 
a surface S enclosed by a circular contour C of radius L » i?ceii as shown in Figure 16. The contribution to the 
flux from cells "inside" the contour is identically zero: the only nonzero contribution comes from cells along C. The 

1/2 

number of such cells is ^ L/Rceii and therefore the flux through C scales as (L/Rceii) ■ Hence the magnetic field, 
averaged over a length scale L is 



Thus, a significant large-scale magnetic field results if the Galaxy is filled with Crab-like regions. 

The proposal that the galactic dynamo is seeded by magnetic fields first created in stars has a potentially fatal 
flaw: Stars form in protostellar gas clouds. Since the presence of angular momentum in a cloud halts its collapse, 
efficient angular momentum transport is necessary if star formation is to occur. Magnetic field can remove angular 
momentum from protostellar clouds, but if the first fields form in stars, then the first clouds will not have any fields. 
Ncvcirtlielcss, stellar field may be important in establishing the galactic dynamo. Bisnovatyi-Kogan, Ruzmaikin, & 
Sunyaev (1973) (see, also Pudritz & Silk 1989) proposed a model in which fields are generated in protostellar clouds 
by a Biermann-type mechanism and dynamo action. These fields facilitate star formation and the stars that form 
then act as sites for rapid field amplification. Finally, supernova explosions and stellar winds disperse stellar fields 
into the ISM thereby seeding the dynamo. 

A common thread between the battery mechanism described in the previous section and fields produced by stars is 
that small-scale fields can provide the seed for a large-scale dynamo. Indeed, the fiuctuation dynamo is very efficient 
at amplifying fields on scales up to that of the largest turbulent eddies {I ~ 100 pc) to equipartition (1 /xG) strengths 
on a relatively short timescale. The projection of this highly tangled field onto the dominant mode of an acu dynamo 
can act as the seed. The amplitude of this component is, in general, a factor of 100 smaller than the rms field thus 
giving a seed field of Bs — 10^* G. Simulations of an auj dynamo seeded by tangled field were carried out by Poedz et 
al. (1993) and Beck, Poedz, Shukurov & Sokoloff (1994). The initial field was chosen to be random on sub-kiloparsec 
scales with an rms strength of 1 /xG. At first, the amplitude of the field declines sharply while its scale length increases 
rapidly. There follows a period of exponential growth to a present day value for the regular field strength of ~ 3 fxC 

4. Active Galactic Nuclei 

Active galactic nuclei (AGN) are promising sites for the production of galactic and extragalactic magnetic fields 
(Hoyle 1969, Recs 1987, Daly & Loeb 1990; Chakrabarti, Rosner, & Vainshtein 1994, Rees 1994). AGN are powered 
by the release of gravitational potential energy as material accretes onto a central compact object, presumably 
a supermassive black hole. There, dynamo processes can amplify magnetic fields on relatively short timescales. 
Moreover, even if the central region of an AGN forms without a magnetic field, one will develop quickly through 
Biermann battery-type mechanisms. Finally, well-coUimated jets can transport magnetic field energy away from the 
central object and into the protogalactic or intergalactic medium. 

The following order of magnitude estimate, due to Hoyle (1969) gives an indication of the field strengths possible in 
an AGN scenario. The rotational energy of a compact object and the material in its immediate vicinity (total mass 
M) can be parametrized as fM(? where / < 1. Equipartition between the magnetic field and the rotational energy 
in the fluid (achieved by differential rotation and/or dynamo action) implies a magnetic field strength 
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where Vc is the volume of the central region. If this field expands adiabatically to fill the volume Vg ~ 100 kpc"^ of 
the Galaxy, a field with strength 
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will result. As an example, Hoyle (1969) considers the values M = 10'' Mq, / = 0.1, Vg = l^f^ cm^, and = 10^ G 
where one finds Bg ~ 10~^ G. 

Our understanding of AGN has improved to the extent that detailed discussions of these objects as potential sources 
of seed magnetic fields is now possible. Daly & Loeb (1990) proposed that the magnetic field in a particular galaxy 
originates during an AGN phase of that galaxy. This scenario presumes that all galaxies (or at least, all galaxies with 
strong magnetic fields) have compact central objects, a conjecture that is now supported by numerous observations 
(see, for example, Magorrian et al. 1998). An AGN phase is characterized by two oppositely-directed, well-coUimated 
jets which transport material into the ISM. Since magnetic fields in the central region are frozen into the jet material, 
they too arc carried into the ISM. 

Daly & Loeb (1990) described, in detail, the interaction of the jet with the ISM (also see Daly 1990). A shock forms 
where this material collides with the ambient gas. At the same time, a blast wave develops propagating perpendicular 
to the jet axis and carrying magnetized material into the so-called shock cocoon. The timcscalc for this process is 
short relative to the lifetime of the galaxy and material from the central AGN engine can reach the outer part of a 
protogalaxy in a time ~ Vg/Lj where Vs is the shock speed and Lj is the length of the jet. An order of magnitude 
estimate for the strength of the galactic field that results is obtained by assuming that the total magncrtic fic^ld energy 
in the jet resides ultimately in the disk galaxy that develops, i.e., VjB^j = VdB^ where Vj and Bj are the volume 
and field strength for the jet and Vd and Bd are the corresponding quantities for the disk. If the jet has a mean 
cross-sectional area ttRj and scale length Lj we have 



where h and Rd are the characteristic height and radius of the disk. Typically, Rj ~h~ 0.5 kpc, Rd — Lj ~ lOkpc 
and Bj ~ 10 /xG. With these values, Daly & Loeb (1990) estimated that the galactic scale would have a strength 
~/iG. 

An essential feature of Daly and Loeb's model is that each galaxy produces its own seed field. An alternative is 
that a population of AGN at high redshift 'contaminates' the protogalactic medium with magnetic field prior to the 
epoch when most galaxies form (Rees 1987, 1994; Furlanetto & Loeb 2001). The jets that originate in AGN often 
end in giant radio lobes which are tens of kiloparsecs in size and have magnetic fields > 10/U G. The field in a typical 
protogalaxy due to pre-galactic AGN is 



where uagn is the number density of AGN. Thus, a relatively small number of AGN can seed the protogalactic 
medium with fields in excess of those produced by battery-type mechanisms. 

The amplification of magnetic fields in an accretion-disk environment has been the subject of a number of studies. 
Pudritz (1981), for example, demonstrated that a mean-field dynamo can operate in a turbulent accretion disk. The 
Balbus-Hawley instability provides an alternate route to strong fields since a weak magnetic field leads to turbulence 
which in turn can rapidly amplify the field (Brandenburg et al. 1995; Hawley, Gammie, & Balbus 1996). 

Kronberg, Lesch, & Hopp (1999), Birk, Wiechen, Lesch, & Kronberg (2000) considered a somewhat different scenario 
in which the seed fields for spiral galaxies such as the Milky Way are created in dwarf galaxies that form at a redshift 
z ~ 10. In their model starburst-driven superwinds rather than jets transport magnetic fiux into the IGM. In the 
hierarchical clustering scenario, dwarf galaxies are the first to form and hence the first to support active star formation. 
In addition, the consensus is that these systems suffer substantial mass loss via supcrnovac and stellar winds. The 
scenario requires that strong magnetic fields exist at early times in dwarf galaxies a reasonable assumption given 
that star-formation and outflow activity increase the amplification rate in models such as Parker's (1992) modified 
aw-dynamo. 

C. Seed Fields from Early Universe Physics 




(111) 




(112) 



The very first magnetic fields may have been created during an early universe phase transition. These events 
typically involve fundamental changes in the nature of particles and fields as well as a significant release of free energy 
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over a relatively short period of time, two conditions that lead naturally to electric currents and hence magnetic fields. 
Often, the question is not whether magnetic fields are created during an early universe phase transition, but whether 
these fields are appropriate seeds for galactic dynamos. 

Causality imposes a fundamental constraint on early universe scenarios for the generation of seed fields. The Hubble 
distance, Lnit) = c/H{t), sets an upper bound on the size of a region that can be influenced by coherent physical 
processes. (In a radiation or matter dominated Universe, Lh is, up to a constant of order miity, equal to the causal 
horizon, i.e., the age of the Universe times the speed of light. Even if an inflationary epoch changes the causal structure 
of the Universe, Lh still sets an upper limit on the scale over which physical processes can operate.) In comoving 
coordinates, the Hubble distance Xh = Lnit) / a{t) reaches the scale \h — 100 kpc of galactic disks at a time t ~ 10''s 
after the Big Bang. Contrast this with the age of the Universe at the electroweak phase transition {Iew — 10~^^s) 
or QCD phase transition {Iqcd — 10~^ s) and the nature of the causality problem becomes clear: Fields generated 
in the early Universe have a coherence length that is much smaller than what is required of seed fields for a galactic 
dynamo. 

Many of the scenarios that operate after inflation rely on either statistical fluctuations of strong small-scale fields 
to yield (weak) large-scale ones or dynamical processes, such as an inverse cascade of magnetic energy, to channel 
field energy from small to large scales. Conversely, scenarios that operate during inflation produce fields on scales up 
to the present-day Hubble radius. However, these models have their own set of difficulties. 

1. Post-inflation Scenarios 

We begin our discussion with the cosmological QCD phase transition (e.g., Boyanovsky 2001). At high temperatures, 
quarks and gluons are weakly coupled and exist as nearly free particles in a plasma. The transition to the hadronic 
phase, in which quarks are bound into mesons and baryons, occurs at a temperature Tqcd — 150 MeV. The order of 
the quark-hadron phase transition is not known. If it is second order, the transformation from quark-gluon plasma to 
hadrons occurs adiabatically, i.e., approximate thermodynamic equilibrium is maintained locally and at each instant 
in time. A more dramatic sequence of events occurs if the phase transition is first order. As the Universe cools below 
Tqcd, bubbles of hadronic phase nucleate and grow. Shocks develop at the bubble walls, latent heat is released, and 
the Universe reheats back to Tqcd- The two phases now coexist with the hadronic regions growing at the expense of 
regions still in the quark phase. 

Hogan (1983) was the first to investigate the possibility that magnetic fields could be generated in a cosmological 
first-order phase transition. His model assumes that battery and dynamo processes create and amplify magnetic 
fields, Bb that are concentrated in the bubble walls. When the walls collide the fields from each bubble are "stitched" 
together by magnetic reconnection. In this way, magnetic field lines, following random paths, can extend to scales 
much larger than the characteristic scale. Lb, of the bubbles at the time of the phase transition. The situation is 
analogous to the one that arises with seed fields from stellar outflows and supernovae. The component of the magnetic 
field that is coherent on scales L ^ Lb has a typical strength — Bb {Lb/ L)^^"^ . 

Detailed calculations of magnetic field generation during the electroweak and QCD phase transitions, under the 
assumption that they are first-order, have been carried out by numerous groups. Quashnock, Loeb, & Spergel 
(1989), for example, demonstrated that a Biermann battery operates during the QCD phase transition. The baryon 
asymmetry of the Universe implies that there are more quarks than anti-quarks. If the number densities of the light 
quarks up, down and strange (charges 2/3, —1/3, and —1/3 respectively) were equal, the quark-gluon plasma would 
be electrically neutral. However, the strange quark is heavier than the other two and therefore less abundant so that 
there is a net positive charge for the quarks. This charge is compensated by an excess of negative charge in the 
lepton sector. Shocks that develop during the nucleation of hadronic bubbles are characterized by strong pressure 
gradients which affect the quarks and leptons differently. Therefore, electric currents develop as bubble walls sweep 
through the quark-gluon plasma. Quashnock, Loeb, & Spergel (1989) estimated the strength of the electric fields 
to be E' ~ (Vpe) /eJi-e ^ eSkT^QY^/LB where 5 = Lb^Pc/Pc characterizes the magnitude of the pressmc gradient 
and e is the fractional difference between energy densities of the quark and lepton fluids. Lb depends on the details 
of the transition. A reasonable estimate is Lb ~ 100 cm which is an order of magnitude or two smaller than the 
Hubble distance at the time of the phase transition. With e ~ (5 ~ 0.1, Quashnock, Loeb, & Spergel found U ~ 5G. 
Lb sets the scale for the coherence length of the field. A distance Lb — 100 cm at tqcD corresponds to a present 
present-day distance I ~ 6 x 10^"^ cm ~ 4AU. The field strength at recombination on these scales is ~ 2 x 10^^^ G. 
Assuming B oc (Hogan 1983) this yields a galactic-scale {L ~ 100 kpc) field with strength at recombination of 

B ~ 6 X 10-32 Q 

At the time of the QCD phase transition, the energy density of the Universe is ^ 1 Gev fm"^ corresponding to 
an equipartition field strength of B^q — 10^^ G. Thus, the battery mechanism proposed by Quashnock, Loeb, & 
Spergel (1989) taps into a tiny fraction of the available energy. Cheng & Olinto (1994) and Sigl, Olinto, & Jedamzik 
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(1997) suggested that stronger field-generating mechanisms operate during the coexistence phase that foUows bubble 
nucleation. As hadronic regions grow, there is a tendency for baryons to concentrate in the quark phase. Essentially, 
the bubble walls act as 'snowplows' sweeping up baryons. In doing so, they create currents of order J ^ en^v where 
n+ is the number density of positive charge carriers in the quark phase and v is the typical velocity of the bubble 
walls. The corresponding magnetic field is B en+vrd where is the thickness of the charge layer. For reasonable 
parameters, this leads to an estimate -Bqcd — 10^ — 10^ G on scales of 100 cm at iqcD- If i? oc L~^/^, a field strength 
at recombination on 100 kpc scales (comoving coordinates) will be i? ~ 6 x 10~^^ — 10~^^ G. 

A first-order electroweak phase transition can also generate large-scale magnetic fields (Baym, Bodeker, & McLerran 
1996; Sigl, Olinto, & Jedamzik 1997). During the electroweak phase transition, the gauge symmetry breaks from the 
electroweak group SU{2)l x U{1)y to the electromagnetism group J7(1)em- The transition appears to be weakly 
first order or second order depending on parameters such as the mass of the Higgs particle (Boyanovsky 2001; Baym, 
Bodeker, & McLerran 1996 and references therein). If it is first order, the plasma supercools below the electroweak 
temperature Tew — 100 GcV. Bubbles of broken symmetry phase nucleate and expand, eventually filling the Universe. 
As in the case of the QCD transition, we write the bubble size Lb = JbLh where /b — lO^'' — 10~^ and Lr — 10 cm. 
The wall velocities are believed to be in the range Uwaii — (0.05 — 0.9)c. Baym, Bodeker, & McLerran 1996 discussed 
the structure of the bubble walls and the associated shocks in detail. The key observation is that the fluid becomes 
turbulent where two walls collide. Fully developed MHD turbulence leads rapidly to equipartition of field energy up 
to the scale of the largest eddies in the fluid, assumed to be comparable to L^. Thus, the field strength on this scale 
is 

B ^ (47re)V^TEw)r|w(^)' 

~ (7 X 10^^ - 2 X 10^^) G (113) 

where e = gt,aT^Y^/2 ~ 4 x lO^^GeVfm"^ is the energy density at the time of the electroweak phase transition. 

Magnetic fields can arise in cosmological phase transitions even if they are second order (Vachaspati 1991). In the 
standard model, electroweak symmetry breaking occurs when the Higgs field <p acquires a vacuum expectation value 
(VEV): {(j)) = Tj. Interactions between <j> and the gauge fields are described by the kinetic energy term in the 
Lagrangian which can be written £kin = D^cjiD'^cj) where = 9^ — igA^j, is the covariant derivative and g is the 
coupling constant. (Gauge group indices have been suppressed.) In principle, and Afj, can conspire to give 

{Dn(f>) = 0. However, since all quantities are uncorrelated over distances greater than the Hubble distance at the 
time of the phase transition, (-D^^), in general, does not vanish (Vachaspati 1991). Dimensional arguments imply 
{D^(j>) ^ rj/S, where ^ is the correlation length for the field at the time of the phase transition. The correlation length 
is actually set by the temperature and is much smaller than the Hubble length: ^ ^ hc/kTEw ^ c/Hew- We are 
interested in the electromagnetic field which is but one component of the complete set of gauge fields A^. Vachaspati 
(1991) finds F^™ ~ g~^rj~'^dfj,4>dv4' which, using dimensional analysis, implies a field strength B ~ 5"^^"^. 

To estimate the field strength on larger scales Vachaspati (1991) assumed that executes a random walk on the 
vacuum manifold with stepsize ^. Over a distance L = where N is a large number, the field will change, on 
average, by N^^^r]~^. Thus, the gradient in the Higgs field is 9(0) ^ rj N^^'^ L^^ = ijN^^^'^S^'^ . The magnetic field 
strength, at the time of the phase transition, will be S ~ BcdN~^ oc L~^. On a galactic scale {L = 100 kpc in 
comoving coordinates; A'' ~ 10^^) at decoupling, we find B ~ 10~^^ G. 

A somewhat different analysis by Enqvist & Olesen (1993) suggests that the magnetic fields decrease with scale 
more slowly than suggested by Vachaspati (1991). They argued that the mean magnetic field satisfies 

(B)=0 (^2)1/2^^ (114) 

where (. . .) denotes an average over regions of size A''^. Thus, in their model, B cx L^^/^. 

Dynamical mechanisms can also lead to an increase in the coherence length of a magnetic field that is produced in an 
early Universe phase transition. Cornwall (1997), Son (1999) and Field & Carroll (2000) have considered the transfer 
of magnetic field energy from small to large scales by a process known as inverse cascade. This process occurs when 
there is substantial magnetic helicity in a fluid (Pouquet, Prisch, & Leorat 1976). Assume that small-scale helicity is 
injected into the fluid in a short period of time. Magnetic energy will shift from small to large scales as the system 
attempts to equilibrate while conserving magnetic helicity and total energy. An inverse cascade in a cosmological 
context leads to the following scaling laws for the field energy and coherence length: 

B^s) = (^^y (^Ay^'B,^^(t.) (115) 
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and 
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where ti is the time when the fields are created. In the above analysis, it is assumed that the inverse cascade operates 
during the radiation-dominated but not matter-dominated phase of the Universe. For illustrative purposes, Field & 
Carroll (2000) considered the evolution of fields created at the electroweak scale. The coherence scale and strength 
of these initial fields arc written as L{ti) = JlcH^^ ~ 0.6/Lcm and i?rms — /sVSTre^ ~ 8 x 10^''' fs G where 
and Jb are dimcnsionlcss parameters. Today, these fields would have a coherence length and scale L(to) = f 3/Lkpc 
and i3ims(io) — 5 x 10~^°/bG. Thus, if the generation mechanism is efficient in producing horizon-sized helical fields 
(/b — fh — 1) then the result will be strong fields on very large scales. However, no compelling mechanism is 
known for generating large net helicity in a cosmic fluid and therefore scenarios based on the inverse cascade must be 
considered highly speculative. 

2. Inflation-Produced Magnetic Fields 

The inflationary Universe paradigm provides both the kinematic and dynamical means of producing a nearly scale- 
free spectrum of energy density perturbations (e.g., Kolb & Turner 1990 and references therein). This feature makes 
inflation an attractive candidate for the production of magnetic fields (Turner & Widrow 1988). 

In most models, inflation is driven by the dynamics of a weakly coupled scalar field known as the inflaton. During 
inflation, the energy density of the Universe is dominated by the vacuum energy of the inflaton (Guth 1981; Linde 1982; 
Albrecht & Stcinhardt 1982). Since this energy density is approximately constant, the Universe is in a nearly dc Sitter 
phase with H{t) ~ constant and a ~ exp {Ht). Exponential growth of the scale factor has two important consequences. 
First, during inflation the physical length L corresponding to a fixed comoving scale A = L{t)/a(t) grows relative to 
the (nearly constant) Hubble distance Lh- By contrast, during the radiation and matter dominated epochs, L grows 
more slowly than Lh- Therefore, a given mode of fixed comoving size (e.g., a single Fourier component of a small- 
amplitude density perturbation) starts out subhorizon-sized and crosses outside a Hubble volume during inflation only 
to reenter the Hubble volume during either radiation or matter dominated epochs (Figure 17). Microphysical processes 
can influence this mode during the flrst subhorizon-sized phase. In particular, de Sitter space quantum-mechanical 
fluctuations continuously excite all massless or very light fields (i.e.. all fields whose Compton wavelength fic/m is 
greater than Lh)- Both the amplitude and coherence length of these fiuctuations are set by the Hubble parameter 
so that the energy density per logarithmic wavenumber bin in these excitations at the time when they are produced 
is de/dlnk ~ HH/L^ ~ TiH^/c'^. Therefore, to the extent that H is constant during inflation (or more generally, a 
power-law function of time) the fiuctuation spectrum is scale-free. 

While all massless fields are excited during inflation, their subsequent evolution varies dramatically depending on 
how they couple to gravity. Fluctuations of the inflaton give rise to an energy-density perturbation spectrum at 
second horizon crossing that is consistent with what is required to explain large-scale structure in the Universe (Guth 
& Pi 1982; Starobinskii 1982; Hawking 1982; Bardeen, Steinhardt, & Turner 1983). In simple inflationary Universe 
models, perturbations have an amplitude at second horizon crossing that is approximately constant with scale. This 
result is consistent with what is called for in hierarchical clustering models and is in agreement with measurements 
of the angular power spectrum of microwave background anisotropics. 

Quantum-produced de Sitter-space fluctuations in minimally coupled scalar fields (i.e., fields that do not couple 
explicitly to the Ricci scalar R) and in gravitons can also be significant (see, Kolb & Turner 1990 and references 
therein). By contrast, since electromagnetic fields are conformally coupled to gravity, the energy density in a fluctu- 
ation decreases as a"** leading to amplitudes that are uninterestingly small. Conformal invariance must therefore be 
broken in order to produce significant primeval magnetic flux. 

In a conformally trivial theory, the field equations are invariant under a rcscaling of lengths at each location in space. 
Likewise, conformally flat spacctimcs, such as de Sitter space and the radiation and matter dominated Robertson- 
Walker models, can be written as (time-dependent) rescalings of Minkowski space. It follows that the form of the 
field equations of a conformally invariant theory in a conformally flat spacetime are time-independent. For example, 
the equations of motion for the magnetic field in a Robertson- Walker space time are: 




(117) 



53 



where 77 is the conformal time related to the clock time t through the expression drj = dt/a(t). For massive fields, 
conformal invariancc is broken by the introduction of a length scale, namely, the Compton wavelength. Conformal 
invariance can also bo broken through couplings to other fields. 

Turner &: Widrow (1988) suggested a number of ways of breaking conformal invariance in electromagnetism: (1) 
introduce the gravitational couplings RA^, B.,j_,^A'^A'^ , RF'^'^F^^, etc; (2) couple the photon to a charged field that 
is not conformally coupled to gravity; (3) c;oiiple the; photon to an axion-like field. Only (1) was considered in any 
great detail. The RA^ terms are the least attractive possibility since they explicitly break gauge invariance by giving 
the photon an effective mass. However, computationally, this case is the easiest to analyze and for a wide range of 
parameters (the coupling constants for the various terms) interesting large-scale magnetic fields can be generated. 
The RF"^ terms are theoretically more palpable but the fields that result are very small. 

Numerous authors have attempted to find more natural and effective ways to break conformal invariance. Ratra 
(1992) calculated the spectrum of magnetic fields produced in a set of infiation models characterized by the inflaton 
potential y((/)) cx exp ((/)). Potentials of this type can be motivated by superstring theory. Ratra (1992) assumed a 
coupling of the infiaton to electromagnetism through a term a exp {(j))F^i,F'^'' and found that fields as large as 10~® G 
could be produced. 

Magnetic fields due to a charged scalar field were considered by Calzetta, Kandus, & Mazzitelli (1998) and Kandus 
et al. (2000). These authors found that charged domains form during infiation which give rise to currents and hence 
magnetic fields during the post-infiation era. A different mechanism has been proposed by Davis et al. (2001) who 
show that the backreaction of the scalar field gives the gauge field an effective mass thus breaking conformal invariance. 
The mechanism is attractive, in part, because it operates in the standard model. Actually, it is the standard-model 
Z-boson that is amplified through its coupling to the Higgs field. As infiation comes to a close, the fiuctuations in the 
Z field are transferred to the hypercharge field (i.e., a linear combination of the Z and photon fields). The Z boson 
acquires a mass at the electrowcak scale leaving behind pure magnetic field. Davis et al. (2001) obtain a magnetic 
field strength of order 10^^^ G on a scale of 100 pc provided certain conditions during reheating are met. 

Garretson, Field, & Carroll (1992) analyzed the amplification of inflation-produced electromagnetic fluctuations by 
their coupling to a PGB. This coupling, which takes the form ~ 0E • B, leads to exponential growth but only for 
modes whose wavelength is smaller than the Hubble radius. No amplification is found for modes outside the horizon. 
The net result is that large-scale magnetic fields of an interesting strength are not produced by this mechanism. 

VI. SUMMARY AND CONCLUSIONS 

It was the late 1940's when the Galactic magnetic field was independently proposed by theorists and detected by 
observers. Since then, galactic and extragalactic magnetic fields have been the subject of intense and fruitful research. 
Nevertheless, fundamental questions concerning their origin, evolution, and nature remain unanswered. 

Magnetic fields have been detected in over one hundred spiral galaxies, in numerous elliptical and irregular galaxies, 
in galaxy clusters, and in the Coma supercluster complex. New instruments such as the planned square kilometer 
array radio telescope will no doubt reveal new magnetic structures. 

It is of interest to note that at present, there is no example of a meaningful null detection of magnetic field in a 
collapsing or virialized system. Conversely, only upper limits exist on the strength of truly cosmological magnetic 
fields. The fact that these limits are several orders of magnitude lower than the strength of galactic and cluster fields 
suggests that magnetic fields arc amplified, if not created, during structure formation and evolution. 

The magnetic fields found in spiral galaxies are unusual in that the strength of the large-scale component is 
comparable to that of the tangled component. By contrast, the fields in ellipticals are random on < 100 pc scales. 
Likewise, cluster fields are tangled on the scale of the cluster itself. The distinction no doubt reflects a key difference 
between spiral galaxies on the one hand and ellipticals and clusters on the other. Namely, the stellar and gaseous disks 
of spiral galaxies arc dynamically "cold" , rotationally supported systems while ellipticals and clusters are dynamically 
"hot" . Evidently, the scale of the largest component of the magnetic fleld in any system is comparable to the scale of 
the largest coherent bulk flows in that system. 

The acij-dynamo is the most widely accepted paradigm for the amplification and maintenance of magnetic fields in 
spiral galaxies. The hypothesis that magnetic fields are continuously regenerated by the combined action of differential 
rotation and helical turbulence is compelling especially in light of the observation that the magnetic structures in 
disk galaxies are in general spiral. One may think of these structures as the MHD analogue of material spiral arms. 
Spiral structure is believed to be a wavelike phenomenon where the crests of the waves are characterized by enhanced 
star-formation activity which, in turn, is triggered by an increase in the local density. Likewise, magnetic spiral arms 
may reflect low-order eigenmodes in a disklike magnetized system. A more direct connection between material and 
magnetic spiral structure is evident in certain galaxies where strong magnetic fields appear in the regions between the 
material arms. The FIR-radio continuum correlation provides further evidence in support of a connection between 
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star formation and large-scale galactic magnetic fields. 

The magnetic fields in ellipticals and clusters require a different explanation. Mergers may play the central role 
in the establishment of these fields since they are likely to be present in merger remnants, typically tidally shredded 
spiral galaxies. The magnetic debris from merger events can act as seeds for subsequent dynamo action. In addition, 
the energy released during a merger event can drive turbulence in the interstellar or intercluster medium. It is unlikely 
that cither ellipticals or clusters will support an aw-dynamo since differential rotation in these systems is relatively 
weak. However, they may support fluctuation dynamos in which turbulence amplifies magnetic fields on scales up to 
the size of the largest eddies in the systems. 

From a theoretical perspective, the greatest challenge in the study of galactic magnetic fields comes from the 
tremendous dynamic range involved. The scale radius and height of a typical galactic disk are of order 10 kpc and 
1 kpc, respectively, while turbulent eddies in the ISM extend in size from subparsec to 100 pc scales. Naive arguments 
suggest that even if initially, magnetic energy is concentrated at large scales, in a turbulent medium, there is a rapid 
cascade of energy to small scales. A mean-field approximation, where velocity and magnetic fields are decomposed 
into large-scale and small-scale components, bypasses this problem. The equation for the large-scale magnetic field, 
known as the dynamo equation, incorporates the effects of the small-scale fields through the a and (3 tensors which, 
in turn, attempt to capture the gross properties of the turbulence (e.g., helicity, spatial anisotropy). 

The dynamo equation for an axisymmetric thin disk can be solved by means of a quasi-separation of variables 
which leads to eigenvalue equations in t, (j), R, and z. The f-eigenvalue gives the growth rate of the magnetic field 
while the 0-cigcnvalue characterizes the symmetry of the field under rotations about the spin-axis of the disk. The 
R equation is similar, in form, to the Schrodingcr equation and its eigenvalue feeds back into into the value of the 
growth rate. The separation-of- variables analysis has yielded a number of encouraging results. Chief among these is a 
demonstration of principle, namely, that disklike systems with a rotation curve similar to that of a spiral galaxy, can 
support a magnetic dynamo. Numerical simulations provide the means to stiidy more realistic models. In particular, 
the effects of a finite disk thickness and deviations from axisymmctry can be explored. These investigations suggest 
ways in which bisymmetric and/or odd party magnetic fields can be excited. 

Unless one is willing to accept magnetic fields as a property of the Big Bang, their existence today implies a violation 
of the MHD approximation at some stage during the history of the Universe. While MHD processes can stretch, twist, 
and amplify magnetic field, by definition, they cannot generate new field where none already exist. Proposals for 
the origin of the first magnetic fields are as varied as they are imaginative. For example, interest in the exotic 
environment of the very early Universe, and in particular cosmological phase transitions, has spawned numerous ideas 
for the creation of seed magnetic fields. A perhaps more appealing set of proposals relies on the ordinary astrophysical 
phenomena that occur during structure formation. Magnetic fields will develop in AGN, stars, and the shocks that 
arise during gravitational collapse. Indeed, rough estimates suggest that AGN and/or an early generation of stars will 
yield fields of strength 10""'^"'^ G on galactic scales. Dynamo action can amplify a field of this strength to microgauss 
levels by a redshift z ~ 2, a result consistent with observations of magnetic fields in high-redshift radio galaxies. 

The astrophysical mechanisms mentioned above were proposed at a time when our understanding of structure 
formation was relatively crude. It is in part for this reason that the creation of seed fields and the dynamo have been 
treated as separate and distinct processes. Indeed, most studies of disk dynamos do not make specific references to 
particular models for seed field production. Likewise, few papers on seed fields follow the resultant fields into the 
dynamo regime. 

Today, semi-analytic models and numerical simulations enable us to study galaxy formation in detail, taking into 

account hierarchical clustering, tidal torques from nearby objects, gasdynamics, and feedback from star formation. 
Moreover, observations of high-redshift supernovae, the CMB angular anisotropy spectrum, and large scale structure 
have pinned down key cosmological parameters such as the densities of baryons and dark matter and the Hubble 
constant. In light of these developments, it may be possible to achieve a more complete description of the origin of 
galactic magnetic fields, one that begins with the production of seed fields and follows smoothly into the dynamo 
regime. 
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FIG. 10 Sequence of events illustrating an axisymmetric odd-parity (AO) aa»-dynamo. (a) Dipole-like poloidal (AO) field, (b) 
Single field line from (a) after it has been stretched by difi^erential rotation, (c) Toroidal component from (b) assuming turbulent 
diffusion and/or rcconncction have decoupled field in upper and lower hemispheres, (d) Plumes created from the toroidal field 
by cyclonic events similar to the one illustrated in Figure 9. (e) Field loops from (d) twisted by the Coriolis effect, (f) Poloidal 
loops created from the those shown in (e). 




FIG. 11 Schematic illustration ol Parker's (1992) galactic dynamo, (a) A magnetic field line is distorted by the Parker 
instability, (b) Loops of magnetic field free themselves through reconnection. 
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FIG. 12 Available time for dynamo action as a function of Qm in spatially flat cosmologies (Oa = 1 — i^m)- Right-hand vertical 
axis gives the amplification factor assuming a growth rate F = 2Gyr~'^. We assume Ho = 70kms"^ Mpc. The four curves m 
each of the plots are for different choices of U, or equivalently, «»: = 25 — solid curve; Zi = W — dotted curve; Zi = 5 — 
dashed curve; Zi = 3 — long-dashed curve. The three plots assume different values for tf. (a) Zf = 0; (b) Zf = 0.4; (c) Zf = 2. 
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FIG. 13 Schematic diagram illustrating the generation of vorticity in a region of crossed pressure and density gradients. The 
open circles represent particles of equal mass separated in y. The gravitational force on the particles, mg is the same. Particle 
1 sits in a higher density region and therefore experiences a smaller pressure gradient force. A velocity field with shear and 
vorticity develops. 
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FIG. 14 Schematic diagram illustrating the Bicrmann battery eflFect. Open cirlcs represent ions and electrons (plus and 
minus symbols, respectively). In (a), physical quantities do not depend on y an electric field develops in the x direction, and 
electrostatic equilibrium is achieved. In (b), the density increases in the y direction and electrostatic equilibrium is no longer 
possible. Instead, this situation leads to a time-dependent magnetic field into the page and a clockwise current as shown. 
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FIG. 15 Experimental set-up showing how the Biermann effect can be produced in the laboratory. Radiation from a laser is 
incident on the target and produces a plasma. This plasma is characterized by density and pressure gradients which are not 
parallel leading to a time-dependent magnetic field. The field produces Faraday rotation in radiation that traverses the plasma 
from left to right. 




FIG. 16 Schematic diagram showing a cross-section of a distribution of magnetic dipole cells in a galactic disk. The circular 
contour C has a radius L and encloses 0[L/ RcmY cells. Only cells on the border contribute to the magnetic flux through C. 
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FIG. 17 Evolution of the physical size for two comoving scales, Ai and A2, and the Hubble radius, c/H in an inflationary 
cosmology. During inflation, the Hubble radius is approximately constant and modes cross outside the horizon. (This event for 
Ai is labelled 'first horizon crossing'. Modes cross back inside the horizon during the post-inflationary epoch. 



